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بهینه یابی سریع خربا با استفاده از الگوریتم ژنتیک موازی تحت قید قابلیت اعتماد 
با بهبود روش جبری نیروها و هوش مصنوعی * 


علی کاوه(؟ وحیدرضا کلات جاری" ۳‏ پوریا منصوریان۳) 


چکیده در بیش‌تر تحقیق‌های انجام‌شده در زمینه‌ی بهینه‌یابی سازه‌ها برای پرهیز از محاسبه‌های طولانی, احتمال خرابی سیستم سازه‌ای 
یم عطیر یحاط یا راز زوی (عض مط تیه هه اس ی ات یاوه از این کرختا اه دزی تب ان یی 
نتایج نامطلوب اقتصادی می‌شود. هدف این مقاله اراثه‌ی چند رامکار برای تسریع بهینه‌یابی حرپا با استفاده از آلکوریتم ژنتیک تحت قید 
قابلیت اعتماد سیسم سازه‌ای اسث, در راهکاز اول در زرش شاخه ورکرائه با موجه به قوائین احنمال, معیاری برای کنا رگذاشتین مداد ی از 
مسیرها و عدم تولید آن‌ها تعیین می‌شود. راهکار دوم بهبود روش جبری نیروها و استفاده از این روش سریع به‌جای روش مرسوم سختی 
است. راهکار سوم استفاده از عامل‌های هوش مصنوعی در شناسایی مسیرهای خحرابی تکراری در تحلیل ه رکروموزوم و هم‌چنین نعیین 
شایستگی به‌ترین کروموزوم‌های هر نسل است که منجر به جلوگیری از محاسبات اضافی و کاهش چشمگیر زمان محاسبات می‌شود. با 
فان زاهکان مور وهای کت ال تضرایی اف ره کین ترجوی دق کر تلو فان محز شم ور یه فان 
شناسایی می‌شوند و تحلی لآن‌ها در همان مرحله و بسیار زودتر از موعد تمام می‌شود. رامکارهای اول, سوم و چهارم با اندکی تغییر قابل 
تعمیر چه قاب‌ها خراهت بر وف آنن حفین به فحیل قابلیت اخماد ر بهیه بان خریلی ۷۵ بل و گنبن ۱۳۰ ميلامع پرواسته فحت گنه در 
نتیجه‌ی استفاده از راهکار اول, در زمان بسیا رکوتاه کران احتمال حرابی باریک‌تری تعیین شد و هم‌چنین استفاده از هر یک از سایر راهکارها 


واژه‌های کلیدی قابلیت اعتماد. احتمال خرابی سیستم سازه‌ای» روش جبری بهبود یافته‌ی نیروها؛ روش شاخه و کرانه . هوش 


مصنوعی. 
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(۱) نویسنده‌ی مسوول: استاده دانشکده عمران دانشگاه علم و صنعت ایران. 
(۲) استادیار دانشکده عمران و معماری, دانشگاه صنعتی شاهرود. 

(۳) استادیار دانشکده عمران و معماری. دانشگاه صنعتی شاهرود. 


مقد مه 

برای اطمینان از عملکرد صحیح یک سازه ابتدا باید 
عواملی را که باعث ایجاد عدم قطعیت می‌شوند. بادقت 
شناسایی کرد. در بسیاری از سازه‌ها بارهای وارد 
جنس مصالح و غیره کمیت‌های غیر قطعی هستند و 
همواره به‌طور تخمینی در مسائل در نظر گرفته 
می‌شوند. ماهیت تصادفی این متغیرها باعث ایجاد یک 
فضای احتمالاتی در روند تحلیل و طراحی سازه‌ها 
قی تو3: 

اولین راه برای در نظر گرفتن این عدم قطعیت‌هاء 
استفاده از ضرایب ایمتی بار و مقاومنت است.: در این 
بجالت: کست‌های,عساله مانید بان رو صاوست: درضتر ایب 
ایمنی نظیر خود ضرب می‌شود و طراحی صورت 
می‌پذیرد. این نوع طراحی اولین سطح طراحی براساس 
نظریه‌ی قابلیت اعتماد سازه‌ها است. این سطح طراحی 
صرفاً برای طراحی یک عضو سازه‌ای منفرد کاربرد 
دارد و اندرکنش اعضای سازه را در نظر نمی گیرد: 

در سطوح پیشرفته‌تر نظریه قابلیت اعتماد (سطوح 
دوم و سوم) کمیت‌های مسأله همراه با واریانس یا تابع 
چگالی احتمال خود در تحلیل و طراحی در نظر گرفته 
می‌شوند. در این حالت ابتدا تعریفی برای خرابی در 
نظر گرفته می‌شود. سپس تابع حالت حدی به‌صورت 
تفاضل مقاومت و بار تعربف می‌گردد. پس از آن 
احتمال خرابی سیستم سازه‌ای با توجه به 
همبستگی‌های میان کمیات مسأله قابل تعیین است. 

بهینه‌یابی سازه‌ها تحت قید قابلیت اعتماد سیستم 
سازه‌ای, بالاترین سطح مسائل (سطح چهارم) در 
زمینه‌ی قابلیت اعتماد است [۱1]. در این حالت هذف 
پیدا کردن سازه‌ای با کم‌ترین وزن است. به‌طوری‌که 
احتمال خرابی سیستم آن از مقدار معینی تجاوز نکند. 

در چند دهه‌ی اخیر محققان متعددی در زمینه‌ی 
قابلیت اعتماد فعالیت داشته و روش‌های گوناگونی 
برای تعیین احتمال خرابی سیستم‌های سازه‌ای ارائه 
کرده‌اند. یکی از متداول‌ترین این راه‌ها؛ تعیین احتمال 


خرابی سیستم سازه‌ای از طریق تعیین مسیرهای 
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بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


خرابی است. روش شاخه و کرانه ( 024 24 طعصعتط 
0 ([2] و روش مسسیرهای کوتساه شده 
(00)عحظ صمتاحتعصصناصه ۳1262060)[ 3] از معروف‌ترین 
این روش‌ها به‌شمار می‌روند. 

در این روش‌ها ابتدا مسیرهای محتمل خرابی از 
مسیرهای کم اهمیت جدا می‌شود و کران بالا و پایین 
احتمال خرابی سیستم سازه‌ای بر اساس مسیرهای 
خرابی محتمل برآورد می‌گردد. استفاده از اين روش‌ها 
با مشکلاتی همراه است که عمده‌ترین معایب آن. زمان 
زیاد لازم برای تعیین مسیرهای خرابی و اختلاف زیاد 
کرانه‌ای بالا و پایین است. برخی از محققان تلاش‌هایی 
برای رفع این مشکلات انجام داده‌اند [4], لیکن هنوز 
امکان استفاده از این روش‌ها در مسائل بهینهیابی 
سازه‌ها با توجه به زمان‌بر بودنشان وجود ندارد. 

بنابراین بسیاری از محققان برای تعیین احتمال 
خرابی سیستم سازه‌ای و استفاده از آن در مسائل 
بهینه‌سازی از فرضیات ساده‌کننده و محافظه‌کارانه 
استفاده کرده‌اند. کریستنسن و موروتسو [2] توغان و 
دالغلو [5,6] و قاسمی. قلعه‌نوی و مستخدمین حسینی 
[۷] با یک فرض ساده‌کننده و محافظه‌کارانه», احتمال 
خرابی کل خرپا را برابر با مجموع احتمال خرابی 
اعضای آن در نظر گرفتند. آن‌ها با استفاده از این فرض 
ساده‌کننده به بهینه‌یابی سازه‌های خرپایی پرداختند. 

بدیهی است استفاده از این فرضیات محافظه کارانه 
در بهینه‌یابی منجر به نتایج نامطلوب اقتصادی می‌شود. 
برای رفع این نقیصه در این تحقیق چهار راهکار ارائه 
کات که با استفاهه اد اوه تصواتب نفک سا 
فرضیاتی منطقی‌تر و واقعی‌تر تعیین می‌شود و از طرف 
فیک تا اه تمو ان رفن گنای 
می‌یابد. 

در راهکار اول برای تعیین احتمال خرابی یک 
مسیر. عضوی که خرابی آن و خرابی عضو اول مسیر 
کوچک‌تسرین احتمال مشترک را پدید می‌آورند» 
شناسایی می‌شود. سپس کلیه مسیرهای بعدی که از 
عضو شناسایی‌شده منشعب می‌شوند. کنار گذاشته 
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می‌شود. عدم تولید این مسیرها باعث صرفه‌جویی 
بسیار قابل توجهی در زمان محاسبات می‌شود. 

در راهکار دوم با بهبود در روابط روش جبری 
نیروها (326000 10۲66 6ن072عع۸1) امکان استفاده از 
آن برای خرپاهای با توپولوژی متغیر ایجاد می‌شود. 
روش جبری نیروها در تحلیل خرپاها در دو مورد نسبت 
به روش مرسوم سختی برتری دارد. اول اين‌که نیروهای 
اعضا به‌طور مستقیم و بدون واسطه تعیین می‌شود و به 
همین دلیل برای تحلیل قابلیت اعتماد کارایی مناسبی 
دارد. مزیت دوم این که روش نیروها نیاز به حل دستگاه 
معادلات با بعد کوچک‌تری دارد (به تعداد مجهولات 
اضافی به‌جای تعداد درجات آزادی فعال در روش 
سختی). بنابراین این روش برای مسائلی مانند بهینه‌یابی 
یر اه قوس تعایل کرت باا یرو ارت 
کارایی بیشتری دارد [9 ,8]. هم‌چنین ماتریس‌های 13 و 
,3 که در این روش برای تعیین نیروهای اعضا به‌کار 
می‌روند (در ادامه تشریح خواهد شد). مستقل از 
خحصوصیات مصالح و سطح مقطع اعضا است. بنابراین 
کافی است فقط یک بار این ماتریس‌ها برای یک خرپا 
محاسبه شوند و پس از آن می‌توان به کرات در تحلیل 
کر تا معط سر ی مورک وی نا تو وان ی 
یکسان از آن‌ها استفاده کرد [9 ,8]. 

برای تعیین احتمال خرابی یک خرپای نامعین 
ساتتلزم باید: یل هاق:مکور حرپا بسن از حتاف شیر 
میله‌ی آن تا رسیدن به ناپایداری صورت پذیرد. 
خرپاهای به وجود آمده در اثر حذف پیاپی میله‌ها از 
لحاظ توپولوژی با یکدیگر تفاوت دارند. 

در این تحقیق از روابط روش نیروهابه گونه‌ای 
استفاده شده است که بتوان از آن در تحلیل خرپابا 
توپولوژی متفاوت استفاده نمود و در نتیجه با 
جلوگیری از تحلیل‌های تکراری. سرعت محاسبات را 
نالا ره 

راهکار سوم استفاده از تکنیک‌های هوش مصنوعی 
است. در این تحقیق در دو فاز متفاوت از عامل‌های 
هوشمند استفاده می‌شود. فاز اول ذخیره‌سازی 
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مسیرهای خرابی و شناسایی مسیرهای تکراری در 
تحلیل قابلیت اعتماد هر کروموزوم است. با شناسایی 
مسیرهای خرابی تکراری و جلوگیری از انجام 
محاسبات اضافی برای آن‌ها؛ در زمان محاسبات 
صرفه‌جویی می‌شود. یکی از محاسن فاز اول راهکار 
موم آیونآمتتی که از هرهاق رای عفن شاه مرای 
کروموزوم اول نسل اول می‌توان بسرای سایر 
کروموزوم‌های نسل‌های بعدی الگوریتم ژنتیک استفاده 
کرد. 

در ف از دوم از عامل هوشمند یادگیرن ده در 
جلوگیری از محاسبات اضافی در تعیین شایستگی 
کروموزوم‌های هر نسل استفاده شده است. در بهینه‌یابی 
استقده از لگوریتم ژنتیک پنن ان دنت چید نسل, 
تعداد کروموزوم‌های تکراری افزایش می‌یابد. حتی 
ممکن است از یک کروموزوم خاص چندین کپی در 
جمعیت یک نسل موجود باشد. با استفاده از عامل 
هوشمند فاز دوم فقط یک‌بار این کروموزوم‌ها تحلیل و 
شایستگی آن‌ها ذخیره می‌شود. و در صورت تکرار این 
کروموزوم‌ها در نسل‌های بعدی. مورد شناسایی قرار 
می گیرند و از تحلیل‌های تکراری جلوگیری می‌شود. 
راهکار سوم در نسل‌های پایانی الگوریتم ژنتیک نقفش 
بسیار موثری دارد. 

در راهکار چهارم از مقدار احتمال خرابی مجاز 
(قید مسأله‌ی بهینه‌یابی) برای کاهش حجم محاسبات 
در تعیین احتمال خرابی کروموزوم‌ها استفاده می شود. 
با استفاده از این راهکار کروموزوم‌هایی که احتمال 
خرابی آن‌ها به مقدار قابل توجهی بیش‌تر یا کم‌تر از 
با متا هوتسن تعانسا رن سای 
آن‌ها در همان مرحله و بسیار زودتر از موعد تمام 
می‌شود. این راهکار در نسل‌های اول الگوریتم ژنتیک 
همع[ هنوز ماهیت تصادفی دارد از کاربرد 
بیش‌تری برخوردار است. 

در ادامه هر یک از راهکارها تشریح می‌شود. پس 
از آن ارتباط بین آن‌ها مورد بحث و بررسی قرار 
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می‌گیرد و در انتها با ذکر مثال‌هایی تأثیر و کارایی آن‌ها 


ارزیابی می‌شود. 


احتمال خرابی سیستم‌های سازه‌ای 

یکی از تعاریف خرابی برای سازه‌ها رسیدن به 
تیدا دای ای اتکی کرت نوا 
یرای کوتا کون ها ا کار شروش لک اسان 
وقوع این مسیرهای خرابی با یکدیگر متفاوت است. 
بسیاری از اين مسیرها احتمال وقوع بسیار کوچکی 
دارد به‌طوری‌که محاسبه‌ی آن‌ها تأثیری در ارزیابی 
اتیال عض ان کر تاره دزد یط وتان فلس | 
افزايش می‌دهد. در نتيجه معمولاً احتمال خرابی 
سیستم‌های سازه‌ای با کنار گذاشتن این مسیرها و در 
نظر گرفتن مسیرهای خرابی محتمل و با استفاده از 
کرانه‌های متداول نظیر کرنل. دیتلفسن یاغیره تعیین 
می‌شود[ 11 ,10]. 

نحوه‌ی تولید مسیرهای خرابی و تعیین احتمال 
آن‌ها با توجه به همبستگی‌های موجود در مراجع 
[2:12] قابل مطالعه است. 
سفق کراکر قز شید فا هریت ام قای اب 
مسیرهای خرابی حاکم و تعیین احتمال خرابی سیستم 
پاهای وی موی اه وکا اس اه فزه 
اشیت: | 2 
عملیات شاخه‌ای. در عملیات شاخه‌ای. عضو (۲) 
در مرحله‌ی ۶ ام از یک مسیر چنان انتخاب می‌شود که 
خرابی آن نسبت به سایر اعضای باقیمانده. بیش‌ترین 
احتمال مشترک را با خرابی اولین عضو آن مسیر یعنی 


1 داشته باشد. 


((۳۵ 6 )هه ِِ (۳ 00 ۳ 


)۱( (وهماحصمجه 1۷۵0نناو) ع [ 


بالا نویس در رابطه‌ی فوق نشان دهنده‌ی 
خرپایی است که اعضای ‏ تا رد خود را از دست 


نشریه مهندسی عمران فردوسی 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


داده است. احتمال وقوع یک مسیر در مرحله ام شامل 
اعضای حذف شده 7 تا ره از رابطه‌ی (۲) به‌دست 


۳ ((۳)]۱۳۶ 2 ۲ 
از 


روند فوق تا مرحله‌ی ۲ که سازه ناپایدار می‌شود 
اقاه ام تیف قفییه یالبرک رای )تاش 
تیان اسف اد سای توا اد 
استفاده می‌شود. در این مقاله از رابطه‌ی (۳) برای 
تنخمین احتمال خرابی هر مسیر استفاده شده است 

[12]. 
دز وی ی ۳ص( .۲ 


۳( ط.ر2-زز (۱۳) 2۳۲۲ 1 


مزیت استفاده از فرمول تقریبی فوق این است که 
نیازی به ذخیره‌سازی داده‌های مربوط به کلیه‌ی مراحل 


به مرحله‌ی اول یعنی سازه‌ی اولیه ذخیره شود. 


تنظیم کران‌های پایین و بالا. بر اساس روش کرنل 
احتمال خرابی سازه‌ی بزرگ‌تر با مساوی حداکثر 
احتمال خرابی کلیه‌ی مسیرها و کوچک‌تر یا مساوی 
مجموع احتمالات خرابی مسیرها است. لذا در روش 
شاخه و کرانه طبق روش کرنل کران بالا و پایین 
احتمال خرابی یعنی و نطبه شرح زیر تعیین 


می‌شود. 


کران پایین. پس از تخمین احتمال خرابی هر مسیر 
جدید ,و,:1در حین عملیات شاخه‌ای (از رابطه‌ی ۳) 
شرط زیر کنترل می‌شود. 
آیا م1 از احتمال محتمل‌ترین مسیری که 
تاکنون به‌دست آمده یعنی ,1۷2۳ بزرگ‌تر است با 
خیر؟ در صورت مثبت بودن پاسخ مقدار ,یم در 
قرار داده می‌شود. 
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63 جوم ۳ < ۷2۳ جح ۷2۳ < تسش 1 
ف ای سا سل ناسا اعسال 
حتف تشرین مین رای اشت یت ب زرابم دز بایسان 
عملیات شاخه و کرانه مقدار ۷12۳7 به‌عنوان کران 
پایین احتمال خرابی سیستم سازه‌ای در نظر گرفته 
می‌شود. 
6( ,۷2۳ 2< ۳۲ 
کران بالا. پس از تعبین بی,.۳ مطابق آن‌چه گفته شد 
این مقدار با مجموع احتمالات خرابی مسیرهایی که 
تاکنون به‌دست آمده یعنی ,50۳ جمع می‌شود. 


0( چیه ۳۰ + مگ ع ,طصت؟ 


همان گونه که بیان شد »برابر با مجموع 
احتمالات مسیرهای خرابی است. بنابراین در پایان 
عملیات شاخه و کرانه ٩۳,‏ به‌عنوان کران بالای 


احتمال خرابی سازه در نظر گرفته می‌شود. 
۹2 ,م٩‏ > ۲۳ 


عملیات کرانه‌ای. این عملیات مشخص می کند که آیا 
اعضایی که قرار است در هر مرحله برای عملیات 
شاخه‌ای انتخاب شوند. منجر به یک مسیر خرابی 
محتمل می‌شوند یا خیر. اعضایی که شرط فوق را 
نداشته باشند توسط عملیات کرانه‌ای کنار گذاشته 
می‌شوند. 

احتمال وقوع یک مسیر در طول عملیات شاخه‌ای 
مرتباً کاهش می‌يابد. به‌عبارت دیگر احتمال آن‌که 
اعضای ۱ تا * خراب شود کم‌تر از احتمال خرابی 
اعضای ۱ تا ۳-1 است. اکنون زمانی که در عملیات 
شاخه‌ای قرار است اعضا در مسیرهای خرابی جدید 
انتخاب شوند می‌توان احتمال وقوع آن‌هارا با 
7 ۵۳0 مقایسه کرد. چنان‌چه احتمال وقوع مسیر 
جدید قبل از تشکیل مکانیزم کوچک‌تر از عدد فوق 
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فاشتته سا ات اقه شین عصی وهای تدش مان آن 
نتیجه گرفت که مسیر پادشده نمی تواند مسیر مناسبی 
عمل کرانه‌ای یا شاخه‌ای می‌گویند. انتخاب مقدار 1 
بستگی به دقت مسأله و نظر طراح دارد. هرچه مقدار 
بزرگ‌تر در نظر گرفته شود تعداد مسیرهای صرف 
نظر شده کم‌تر بوده و دقت محاسبات بالاتر. از طرف 
دیگر می‌توان نتیجه گرفت که محتمل ترین مسیری که 
در عملیات کرانه‌ای کنار گذاشته شده است. احتمال 
وقوع کوجک‌تری نسبت به ۲ ۷۵۳,۰10 دارد. 
عملیات شاخه و کرانه تا جایی ادامه می‌یابد که دیگر 


هیچ عضوی برای انتخاب وجود نداشته باشد. 


روش جبری نیروها در تحلیل استاتیکی سازه‌ها 


روابط اصلی و تاریخچه‌ی روش جبری نیروها را 
می‌توان در مراجع [13,14] پیدا نمود. کاوه و کلات 
جاری در سال‌های ۲۰۰۲ و ۲۰۰۶ در مراجع [8,9] با 
ارائه‌ی فرمول‌بندی به‌تر در روش جبری نیروهاء از آن 
در بهینه یابی خرپاها توسط الگوریتم ژنتیک استفاده 
کردند. پس از آن رهامی, کاوه و قلی‌پور [15] از روش 
مجتمع نیروها و روش‌های انرژی در بهین»‌سازی 
توپولوژی خرپاها توسط الگوریتم ژنتیک استفاده 
کردند. 

در روش نیروها معادلات تعادل کل سازه. مستقل 
از خصوصیات مکانیکی مصالح به‌صورت رابطه‌ی 
ماتریسی زیر قابل بیان است : 


۳ (۸) 


در رابطه‌ی بالا 2 برابر درجات آزادی فعال گرهی, 
0 برابر تعداد اعضاء ۳ بردار 9 بعدی بارهای گرهی, ۳ 
ماتریس تعادل ۳ و ۲بردار 7 بعدی نیروهای 
عضوی است. هر یک از سطرهای ماتریس 11 رابطه‌ی 
تعادل نیروهای اعضای متصل به گره مربوط با بارهای 


خارجی در جهت درجه‌ی آزادی مفروض را بیان 


نشریه مهندسی عمران فردوسی 


می‌کند. به اين ترتیب برای یک گره فعال خرپا در 
فضای سه‌بعدی. سه شرط 20<,ر(. 20 ,؟ر2 و 
0 ,2۲ ارتباط نیروهای اعضای متصل به گره مورد 
نظر با بارهای خارجی را برقرار نموده‌است و سه سطر 
پیایی از ماتریس 11 را پر می‌کند. 

در یک سازه‌ی پایدار صسلب بهازای 
20 1-7-7 ماتریس 3 دارای مرتبهی کامل 
سطری است؛ به عبارت دیگر رنک 3 برابر با 9 است. 
1 تعداد مجهولات اضافی يا درجات نامعینی استاتیکی 
سازه می‌باشد. 

روش نیروها مشتمل بر دو گام در نظر گرفته 
می‌شود. در گام نخست ماتریس‌های ,13 و ,ظ به‌نحوی 
پیدا می‌شود که بردار نیروهمای عضوی ۲ در رابطه‌ی 
(۸) به‌صورت زیر بیان شود: 
0( ۳" 


او 


در این رابطه مظ یک ماتریس 10*8 بعدی است 
به‌طوری‌که ۲۲۰ ,3 معادل ماتریس واحد 8 
می‌شود. هر ستون مق[ مبین نیروهای ایجاد شده در 
اعضای سازه‌ی صلب معین اولیه در صورت اعمال بار 
واحد در جهت هر یک از درجات آزادی گرهی است. 
,3 ماتریس 2*۲ از مرتبه‌ی کامل ستونی است 
به‌نحوی که 13۰ ,3 معادل ماتریس صفر می گردد. تعداد 
پکرهای باس اخشرانی روعات انش ز هایگ 
است و هر ستون آن نشانگر سیستم نیروهای عضوی 
است که با خود (در غیاب بارهای خارجی) در حال 
تعادل است و به بردار پوچی با سیستم خحودمتعادل 
معروف است. در رابطه‌ی )٩(‏ بردار و از مرتبه‌ی ] به 
بردار نیروهای مجهول اضافی و ماتریس 1 به ماتریس 
پایه‌های پوچی با ماتریس تنش خودمتعادل موسوم 
است. 

جهت تعیین روابط کلی برای ماتریس‌های ,13 و 
ماتریس تعادل ]7 رابطه‌ی (۸) به دو زیرماتریس 
,1 و را تجزیه می‌شود. در یک سازه‌ی صلب. 
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زیرماتریس ,3 دارای مرتبه‌ی کامل سطری 0 است (به 
عبارت دیگر ",11 وجود دارد) و مربوط به نیروهای 
داخلی اعضای سازه‌ی معین اولیه صلب یعنی و« 
می‌باشد. زیرماتریس 11 نیز دارای مرتبه کامل ستونی ] 
بوده و اختصاص به نیروهای مجهول اضافی م6 
دارد. به‌این ترتیب رابطه‌ی (۸) به‌صورت جدید زیر 
قابل نمایش است: 


۰ رات | و ]۴ 


چنان‌چه رابطه‌ی فوق برای بردار نیروی ۳ حل 
شود نتیجه می‌شود: 


0 وم ]+۲۰۳ 


حال اگر به رابطه‌ی فوق بردار 4 افزوده گردد. در 
قیاس با رابطه‌ی )٩(‏ می‌توان روابط زیر را برای 
ماتریس‌های وق و رظ نتیجه گرفت: 


9 ۳ ن 2 و رن 


با توجه به این که درایه‌های ماتریس‌های ۳ و و با 
توجه به شماره‌گذاری اعضا در سطرهای مربوط در 
بردار نیروهای داخحلی ۲ واقع شده‌اند. لذا سطرهای 
ماتریس صفر ,0 در تعریف ,13 به‌طور نظیر 
سطرهای ماتریس واحد ,1 در تعریف 13 بایستی 
به ترتیب در سطرهای مربوط به اعضای مجهول اضافه 
(واقع شده در قرار گیرند. این جابه‌جایی سطرها در 
بعضی نمادهای فرمول‌بندی روش یروها توسط 
ماتریس‌های تعا رکشت نمایش داده شده‌اند. باید توجه 
داشت که بسته به انتخاب مجهولات اضافه در بردار 0 
تعداد متعددی ماتریس ,3 و ,3 می‌توان تشکیل داد. 
روش‌های متنوعی وجود دارد که بتوان ماتریس‌های ,13 


ور خوب که پر صفر نواری و موزون هستند. 
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به‌دست آورد [13]. 
در گام دوم روش نیروها با استفاده از اصل کار 
مجازی و سازگاری تغییر شکل‌ها؛ دستگاه معادلات زیر 


برای کسب بردار 4 حل می‌شود : 


۱۳( 5,5,۲ )-- 3۳,5۰ ) 
در رابطه‌ی فوق ,] ماتریس قطری 10*1 
مربوط به نرمی اعضا و 1 ,۳ 13<) ماتریس نرمی 
کل سازه می‌باشد. یک ماتریس 13 نواری سبب نواری 
شدن ماتریس 0 خواهد شد [13]. تغییر مکان‌های 
گرهی نیز توسط رابطه‌ی زیر قابل ارزیابی است: 


13 مه 18 


یر رو 0 


در این تحفیق از روش لآ برای تعیسین 
ماتریس‌های و و ر,ظ استفاده شده است. در این روش 
ماتریس نرمی ,13 ,۲ 13 < 6 متقارن و معین- مثبت 
به‌دست می آید و با استفاده از روش چولسکی می‌توان 
بردار 4 را به‌سرعت تعیین نمود و نیازی به کنترل‌های 
وجود ندارد. روش تجزیه گاوس- جردن 
(۱۵۵0۵ جمتازومم۵۵ع۲(۵ مقل10۳ - ووناع)) و روش 
باز گشتی (۳۵0۱0۲12۵100 نان »261 صت؟ تان) نیز در 


مرجع [۱۱۹ معرفی نله اننتست. 


روش افراز لاب 
فرمول‌بندی اولیه‌ی این روش توسط کانکو و همکاران 
[17] صورت پذیرفت و پس از آن کاوه و کلات جاری 
در مرجع [8] با ارائه‌ی روابط جدید. راندمان و سرعت 
ان را افزایش دادند. در این روش ماتریس تعادل 
ب‌صورت [یا,11 ]11 به‌نحوی تجزیه می‌شود که 
ند ی ی و موی 
ستون‌های آن منتخبی از ستون‌های غیروابسته (مستفل) 
ماتریس 1 می‌باشد. م1 ماتریس غیرویژه‌ی پایین مثلشی و 
ماتریس غیرویژه‌ی بالا مثلئی با درایه‌های واحد بر 
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قطر اصلی آن می‌باشد. با استفاده از روابط (۱۵) 
ستون‌های ماتریس‌های بآ و تا به‌طور هم‌زمان و یکی 
پس از دیگری تولید و در صورتی که درایه‌ی قطر 
اصلی ماتریس ,آ صفر شود برای جلوگیری از ویزه 
شدن ماتریس را ستون نظیر از ماتریس 1 که یک 
ستون وابسته است و متناظر با شماره‌ی عضو مجهول 
اضافی. نادیده گرفته می‌شود و محاسبات از ستون 
بعدی دنبال خواهد شد. به‌این‌ترتیب» ستون‌های وابسته 
از ماتریس ]7 که یکسان با تعداد مجهولات اضافی 


به‌ازای ...وش و[ <-[ٌ 
1۳-1 
انس م2 و 
1- .ور 2ب < 0 : وا 
1- 
۱ < 1 تلایا 2 - ط ‌ ۳ 
1<ع1 


در روابط بالا [نط [نا و [ننا به ترتیب درایه‌های 
سطر 1و ستون ژاز ماتریس‌های 11 او لا می‌باشند. 
با تشکیل ماتریس‌های تالا ,11 و بل در یک 
سازه‌ی پایدار طبق روابط (۱۲) ماتریس‌های 3 و 3 
به‌دست ی | تلا 
ثِ"- 2 
1 


(۱۲ تکراری) 8 


",را می‌توان به‌راحتی با استفاده از روش‌های 
محاسبات عددی مانند روش تجزیه ناوید 
[168]. زیرماتریس ۵ در ماتریس 3 صفر از مرتبه‌ی 
ای ی ما ی وان از 
مرتبه‌ی ۴( است. سطرهای این دو زیرماتریس پس 
از انجام عملیات جابه‌جایی سطری مناسب. به‌ترتیب در 
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سطرهای مربوط به مجهولات اضافی قرار می‌گیرد. 
ماتریس 1 که به‌این‌ترتیب حاصل می‌شود. به‌صورت 
پلکانی است که درایه‌های صفر در گوشه‌ی سمت چپ 
و پایین ماتریس قرار می‌گیرد. مرحله‌ی بعدی 
محاسبات. به کارگیری روابط (۱۳و )٩‏ برای محاسبه‌ی 
بردار مجهولات اضافی و نیروهای داخلی اعضا است. 
تغییر مکان‌های گرهی را می‌توان از رابطه‌ی (۱۶) 
محاسبه کرد لیکن در تحلیل قابلیت اعتماد کاربردی 


ندارد. 


الگوریتم ژنتیک موازی 

الگوریتم ژنتیک ساده توسط مرجع [19] پایه‌ریزی شد 
و طی سال‌های بعد تحقیقات بسیار زیادی درباره‌ی آن 
انجامررگرفت [20-23]. این روش بهینه یابی علی‌رغم 
تمامی مزایایش, بسیار وقت گیر است. خوشبختانه 
روابط حاکم در الگوریتم ژنتیک به گونه‌ای است که 
می‌توان از پردازش موازی در آن استفاده نمود. 

تاکنون تحقیقات زیادی در زمینه‌ی موازی کردن 
پردازش در الگوریتم ژنتیک صورت پذیرفته است. 
پردازش موازی در اکثر این تحقیقات به دو صورت 
عمده انجام شده است [24-26]: موازی‌سازی مستقیم 
[27] و موازی‌سازی جزیره‌ای [28,29]. 

سیستم اول که به ۱۷۵5007-81276 موسوم است. 
شامل یک پردازشگر اصلی و ۰-1 پردازشگر فرعی 
ست. در این سیستم پردازشگر اصلی کروموزوم‌های 
یک نسل را بین پردازشگرهای فرعی تقسیم می‌کند. 
ین پردازشگرها به‌طور موازی به محاسبه‌ی شایستگی 
کروموزوم‌های خود پرداخته و نتایج را به پردازشگر 
صلی ارسال می‌کنند. اکنون وظیفه‌ی تولید نسل جدید 
پر عهده‌ی پردازشگر اصلی است. در اغلب مسائل 
مهندسی حجم بالایی از محاسبات (بین ۸۵ تا ۹۵ 
درصد) صرف تعیین شایستگی هر کروموزوم می‌شود 
1 لذا استفاده از پردازش‌گرهای مسوازی تتاتیر 


قابل توجهی در کاهش زمان محاسبات دارد. 
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بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


در سیستم جزیره‌ای که به الگوریتم ژنتیک 
توزیع‌شده معروف است. ابتدا جمعیت اولیه به چند 
زیرجمعیت (جزیره) تقسیم می‌شود و هر یک از جزایر 
توسط یک پردازشگر به‌طور مستقل به بهینه‌سازی 
می‌پردازد. پس از گذشت چند نسل برای جلوگیری از 
ایجاد بهینه‌ی محلی و ایجاد تنوع در جواب به‌ترین 
کروموزوم‌های هر جزیره به جزایر دیگر مهاجرت 
می‌کنند. در هنگام مهاجرت این کروموزوم‌ها. یک کپی 
از آن‌ها در جزیره‌ی شتا بت که رن سوق 
کروموزوم‌های مهاجر که شایسته‌ترین کروموزوم‌های 
جزیره‌ی مبداً بوده‌اند. جایگزین بدترین کروموزوم‌های 
جزیره‌ی مقصد می‌شوند. عمل مهاجرت به‌صورت 
متناوب تا رسیدن به شرط پایان بهینه‌یابی انجام 
ی تو9: 

به تعداد نسل‌های بین دو مهاجرت. وقفه‌ی 
مهاجرت و به تعداد کروموزوم‌های مهاج نرخ 
مهاجرت می‌گویند. 

عیب عمله‌ی سیستم 51276 -۷125]0 نسبت به 
پردازش جزیره‌ای» حجم بالای داده‌های انتقالی بین 
پرشا کت کر هام فا عیدوت اس ای 001 ]هن اس 
سیستم در هر نسل باید تمامی کروموزوم‌هابه 
پردازشگرهای فرعی ارسال و پس از تعیین شایستگی 
پس فرستاده شود. در حالی که در سیستم پردازش 
جزیره‌ای عمل مهاجرت بین جزایر با تعداد محدودی 
انجام می‌شود. به‌طور کلی پردازش موازی جزیره‌ای 
دارای مزایای زیر است : 

حجم داده‌های انتقالی با توجه به محدود بودن 
تعداد مهاجرت‌های بین جزایر پایین است. 

با توجه به کوچک بودن جمعیت هر جزیره 
عملیاتی مانند مرتب‌سازی کروموزوم‌ها؛ انتخاب و غیره 
سریع‌تر انجام می‌شود. 

مهاجرت کروموزوم‌های شایسته‌ی هر جزیره به 
جزایر دیگر باعث می‌شود تا جواب‌های متنوعی ابجاد 
شود و از ایجاد بهینه‌ی محلی خودداری گردد. 
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در این سیستم امکان استفاده از پارامترهاو 
روش‌های گوناگون الگوریتم ژنتیک در هر یک از 
جزایر وجود دارد. این پارامترها تأثیر قابل‌توجهی در 
بهینه‌یابی دارند. در یک الگوریتم ژنتیک ساده یا یک 
سیستم ۷/25067-512۷6[ لازم است از پیش حندین باو 
بهینه یابی انجام شود تا مقداری صحیح برای این 
پارامترها و روشی مناسب تعیین شود. در حالی که در 
روش جزیره‌ای با توجه به تعدد جزایر امکان 
به کارگیری پارامترها و روش‌های گوناگون وجود دارد. 
لذا تعداد مراحل آزمون و خطا و احتمال گیر افتادن در 
بهینه‌ی محلی به طور قابل توجهی کاهش می‌یابد 
[ 28,29,31]. 

به‌طور خلاصه. هر دو سیستم پردازش موازی 
مستقیم و جزیره‌ای باعث افزایش سرعت محاسبات 
می‌شوند. لیکن روش جزیره‌ای علاوه بر آن که سریح‌تر 
اسیته سرت همگرآیی بیرق داشسته و ثسر اخسال 
قرار گرفتن در بهینه‌ی محلی را کاهش می‌دهد. 

در مرجع [29 ,28] از روش موازی جزیره‌ای برای 
بهینه پایی خریاها استفاده شده است. در تحقیق یادشده 
مقادیر مختلفی برای نرخ ادغام و نرخ جهش در جزایر 
استفاده شده است. در مرجع [۳۱] برای هر یک از 
جزایر پارامترها و قوانین به‌دست آمده از تحقیقات 
معروف گذشته مورد استفاده قرار گرفته است. 

در این تحقیق از روش موازی جزیره‌ای برای 
بهینه‌یابی استفاده می‌شود. بدین منظور یک شبکه‌ی 
کامپیوتری شامل ۵ کامپیوتر با مشخصات جدول (۱) 
به‌کار گرفته می‌شود. هر کامپیوتر به یک جزیره 
اختصاص داده می‌شود. هر یک از جزایر با پارامترها و 
روابط مختص خود به بهینه‌یابی می‌پردازد. اطلاعات 
مربوط به هر جزیره و پارامترهای مربوط در جدول 


( او اه یهاگ 
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جدول ۱ مشخصات سیستم‌های کامپیوتری مورد استفاده 
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در اکثر تحقیقات انجام شده از دو روش برای انتخاب 
جزیره‌ی مقصد استفاده شده است [130. در روش اول 
جزیره‌ی مقصد به‌صورت رندوم انتخاب می‌شود. در 
روش دوم مهاجرت کروموزوم‌ها از هر جزیره به 
جزیره‌ی همسایه صورت می‌پذيرد. ایراد این دو روش 
ان تیه که عیکتخ تبرت از غزایر افتساه سنا 
یکدیگر ارتباط برقرار نکنند. در این تحقیق برای رفع 
ابش تهب الکتوریم یشک (۱) سرا 
مهاجرت کروموزوم‌ها پيشنهاد می‌شود. 


شکل ۱ چگونگی انتخاب جزیره‌ی مقصد 
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در شکل (۱) حلقه‌های درونی و بیرونی به‌ترتیب 
نشان دهنده‌ی جزایر مبداً و مقصد در اولین مهاجرت 
هستند. مطابق این شکل در مهاجرت اول 
کروموزوم‌های شایسته‌ی جزیره 1 ام به جزیره‌ی 1۳1 ام 
مهاجرت می‌کنند. در هنگام مهاجرت دوم حلقه‌ی 
بیرونی در خلاف جهت حرکت عقربه‌های ساعت 
می چرخد و جزیره‌ی 1۳2 ام به‌عنوان جزیره‌ی مقصد 
کروموزوم‌های شایسته جزیره‌ی 1 ام انتخاب می‌شود. 
در 7 امین مهاجرت جزیره 1 ام روبه‌روی خود قرار 
می‌گیرد. انتقال از ییک جزیره به خود آن جزیره 
مفهومی ندارد نذا حلقه‌ی بیرونی یک بار دیگر 
می‌چرخد و جزیره 1*1 ام انتضاب می‌شود. 
به‌این ترتیب» حلقه‌ی بیرونی یک دور کامل می چرخد 
و همه چیز از ابتدا آغاز می‌شود. استفاده از اینن 
الگوریتم باعث می‌شود تا انتقال بین جزایر به صورت 
مطمئن‌تر و سریع‌تری صورت پذیرد. 

وقفه‌ی مهاجرت و نرخ مهاجرت تأثیر مهمی در 
یافتن بهینه‌ی کلی دارد. چنان‌چه عمل مهاجرت زیاد 
صورت پذیرد و تعداد کروموزوم‌های مهاجر زیاد 
باشد. جزایر فرصت تولید کروموزوم‌های تکامل يافته‌تر 
را از دست می‌دهند. در این حالت کروموزوم‌ها 
پیوسته از بین می‌روند و کروموزوم‌های جزایر دیگر 
جای آن‌ها را می‌گیرند. چنان‌چه وقفه‌ی مهاجرت و 
نرخ آن کوچک باشد جزایر نسبت به یکدیگر ایزوله 
می‌شوند و مزایای استفاده از یک جمعیت بزرگ از بین 
می‌رود [25]. 

راهکارهای پیشنهادی 

در این مقاله جهار راهکار برای کاهش حجم و زمان 
محاسبات در فرآیندربهینه‌یابی سازه‌های خرپایی با 
استفاده از الگوریتم ژنتیک تحت قید قابلیت اعتماد 
ارائه شده است: 
- راهکار پیشنهادی در قابلیت اعتماد. 
- بهبود فرمول‌بندی روش جبری نیروها. 
- استفاده از تکنیک‌های هوش مصنوعی. 
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- راهکار پیشنهادی در بهینه‌یابی. 


راهکار پیشنهادی در قابلیت اعتماد 
همان گونه که بیان شد. در این تحقیق احتمال یک مسیر 
خرابی مانند اباب بآ طبق رابطه‌ی تقریبی ۳ برابر 
با کوچک‌ترین احتمال مشترک خرابی اعضای [. و 1 


خط پر نمایش داده شده است. 


شکل ۲ طرح شماتیک از مسیرهای خرابی منشعب از [س[ 


بنابراین اگر کوچک‌ترین احتمال مشترک این 
اعضا مطابق نمودار ون شکل (۳) مربوط به اعضای 1و 
ز باشد آن‌گاه احتمال مسیر خرابی اباب م1 یعنی 
(,۱۳,۲۱۳,/۱۳] ۳0۳ با احتمال مشترک خرابی اعضای 
و ژیعنی (۲۱5 )۳ تقریب زده می‌شود. اين تقریب 
محافظه‌کارانه است زیرا طبق شکل (۳) رویداد 
زیر مجموعهی رویداد ۲15 


است. 


شکل ۳ نمایش کوچک‌ترین اشتراک در نمودار ون 
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این مسأله را می‌توان به‌راحتی به کلیه‌ی شاخه‌های 
یت اشنا از مس تعمیم داد لذا همه‌ی این مسیرها 
زیرمجموعه‌ی رویداد ۳۲۱۳ هستند. این شاخه‌ها با 
نقطه چین در شکل (۲) نمایش داده شده است. 

اکنون اگر طبق رابطه‌ی (۳) رویداد مسیر خرابی 
اما زب با رویداد ,۱۲ 7 تقریب زده شود. آن‌گاه 
احتمال اين رویداد یعنی (۳)۳:117 طبق کران بالای 
کرنل با مجموع احتمال مسیرهای قبلی جمع زده 
ی کی 3: 

در این حالت باید از تولید مسیرهای خرابی دیگر 
که از سا منشعب می‌شود. جلوگیری کرد. زیرا 
همان گونه که گفته شد» همه‌ی آن‌ها زیرمجموعه‌ی 
رویداد ,۱۳ هستند و اجتماع چند رویداد که همگی 
زیرمجموعه‌ی یک رویداد هستند نیز زیر مجموعه آن 
رویداد است. تولید این گونه مسیرها نه تنها چیزی به 
فضای خرابی کل اضافه نمی‌کند بلکه باعث می‌شود 
احتمال آن‌ها بر اساس کران بالای کرنل به‌صورت 
اشتباه با هم جمع شود و مقدار بسیار محافظه کارانه‌ای 
برای کران بالای احتمال خرابی سیستم سازه‌ای تعیین 
گردد. علاوه بر آن تولید این گونه مسیرها حجم 
محاسبات را بی‌دلیل افزایش می‌دهد. 

با حذف مسیرهای منشعب باقیمانده نه تنها کران 
فوقانی کوچک‌تری برای احتمال خرابی سیستم سازه‌ای 
تعبین می‌شود. بلکه سرعت محاسبات نیز افزایش 
می‌یابد. 

با جایگزینی رویداد خرابی یک عضو خرپا با 
رویداد ایجاد یک مفصل پلاستیک. راهکار پیشنهادی 
قابل تعمیم به سازه‌های قابی است. 


بهبود فرمول‌بندی روش جبری نیروها 
در روند تشکیل یک مسیر خرابی زمانی که عضوی از 
خرپا حذف می‌شود علاوه بر توزیع مجدد نیروها باید 
ماتریس سختی کل خرپا اصلاح شود. این امر با صفر 
در نظر گرفتن سختی عضو مورد نظر و تصحیح 
دوآنه‌هایی ان ماتریی خی کت کته چا آن عفر دز 
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جدید تحلیل قابلیت اعتماد شود تا احتمال خرابی اعضا 
ت‌دست ایک و ان آن‌هنا عضو ی را که سین مربرخ 
احتمال خرابی دارد انتخاب شود. این عملیات حذف 
اعضا و تحلیل سازه‌های جدید باید تا وقوع خرابی یا 
همان صفر شدن دترمینان ماتریس سختی کل خرپا 
ادامه یابد. 


1 


-<0 (۱۹ 


گ 

بالائویس ۲ نشان‌دهنده‌ی کاهش به‌وجود آمده در 
ماتریس سختی کل است. هر بار تحلیل قابلیت اعتماد 
خرپا مستلزم حل مکرر دستگاه معادلاتی با تعداد 
درجات آزادی فعال گرهی است. 

در روش متعارف جبری نیروها نیز باید مشابه 
آن‌چه در روش سختی بیان شد. ماتریس 11 را تصحیح 
کرد. پس از آن همه‌ی عملیات از جمله افراز ماتریس 
7 با روش تجزیه لاب تعیین ماتریس‌های 130 و ,13 و 
حل رابطه‌ی (۱۳) را دوباره باید انجام داد. عملیات 
فوق زمان زیادی صرف می کند به‌طوری‌که استفاده از 
روش نرمی منجر به صرف زمان بیش‌تری نسبت به 
روش سختی می‌شود. 

اک نمی توا اعمال تق ای تفن توس تا 
جبری نیروها زمان لازم برای تحلیل خرپا را چنان 
کاهش داد که سریع‌تر از روش سختی به جواب برسد. 

نکته‌ی اصلی این راهکار در این‌جا است که 
ماتریس‌های 30 و 3 ابدا ارتباطی به سختی یانرمی 
اعضا ندارند و در صورت حذف یک عضو کافی است 
درایه‌ی نظیر آن در ماتریس نرمی ,۳ اصلاح شود. با 
اصلاح ماتریس یاد شده کافی است یکبار دیگر 
رابطه‌ی (۱۳) حل شود تا مقادیر نیروهای اعضا تعیین 
شود. به بیان دیگر ماتریس‌های یادشده فقط یکبار 
برای سازه‌ی دست نخورده‌ی کروموزوم اول نسل اول 
تشکیل می شوند و تا انتها بدون تغییر باقی می‌مانند. 

روش بهبود یافته‌ی نیروها دو برتری نسبت به 


روش سختی دارد. اول این که در این روش حل 
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دستگاه معادلات خطی از تعداد درجات آزادی فعال 
گرهی به تعداد درجات نامعینی استاتیکی (که معمولا 
کشت ) تقلیان م‌بانگا دبا توس به این کهساژلات 
فوق در تحلیل قابلیت اعتماد سازه‌ها با دفعات بسیار 
یات سل میسرت انتاده او روشن پرشتهآدی عی تران 
باعث صرفه‌جویی چشمگیر در حجم محاسبات شود. 
مزیت دوم این است که نیروهای اعضا به‌طور مستقیم 
و قبل از محاسبه‌ی بردار تغییر مکان‌های گرهی تعیین 
می‌شود و باعث کاهش بیش‌تر مان محاسبات خواهند 
شد. 

دو نکته در اینجا وجود دارد. اول این که درایه‌های 
ماتریس نرمی را نرمی اعضای خرپا تشکیل می‌دهد. 
نرمی عضو ۶ ام خرپا از رابطه‌ی زیر به‌دست می‌آید. 

۳ 


۱ 2 
۳ 


منجر به بی نهایت شدن نرمی آن می‌شود. در عمل در 
صورت حذف هر عضو از یک عدد بسیار بزرگ برای 
نرمی آن استفاده می‌شود. 

نکته‌ی دوم این که در تحلیل خرپاهای بزرگ زمان 
زیادی صرف انجام ضرب‌های ماتریسی موجود در 
رابطه‌ی (۱۳) می‌شود. در این حالت مشابه تعدیل 
بیان‌شده در روش سختی. می‌توان کنترل کرد که تخیر 
نرمی عضو ام بر چه درایه‌هایی از ماتریس‌های 
حاصل ضرب رابطه‌ی (۱۳) تأثیر دارد. این امر می‌تواند 
صرفه‌جویی قابل توجهی در زمان محاسبات ایجاد کند. 

شایان ذکر است که در روش پیشنهادی. ایجاد 
نیرو در هر یک از میله‌هایی که حذف شده است 


به‌منزله‌ی ناپایداری خرپا خواهد بود. 


استفاده از تکنیک‌های هوش مصنوعی 
در این مقاله از عامل‌های هوشمند در دو فاز متفاوت 


استفاده می‌شود. فاز اول شامل ذخیره‌سازی مسیرهای 


نشریه مهندسی عمران فردوسی 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


خرابی هر کروموزوم و شناسایی این مسیرها است. 
عامل هوشمند در صورت تکرار این مسیرهاء نه تنها در 
ادامه‌ی تحلیل آن کروموزوم. بلکه در تحلیل همه‌ی 
کروموزوم‌های نسل‌های بعدی آن‌ها را شناسایی 
می‌کند. حسن شناسایی این مسیرهای خرابی» همان‌طور 
که از نام آن‌ها پیدا است» این است که شرط خرابی 
قطعاً در مورد آن‌ها برقرار است و نیازی به کترل 
نیست. بنابراین زمان لازم برای کنترل شرط خرابی در 
مورد تمامی این مسیرهای تکراری ذخیره می‌شود. 

در روند بهینهیابی باالگوریتم ژنتیک 
کروموزوم‌هایی که شایستگی بیش‌تری دارند. به کرات 
تولید می‌شوند. عامل هوشمند فاز دوم شایستگی این 
گونه کروموزوم‌ها را ذخیره می‌سازد و پس از تکرار 
آن‌ها مقدار شایستگی ذخیره شده را اعلام می‌کند. به 
این ترتیب از تحلیل اضافی جلوگیری می گردد و حجم 
لها از مخاسیات کار کلاشته مشود 

در ادامه نحوه‌ی کارکرد این دو عامل به‌طور 


عامل هوشمند فاز اول 

در یک خرپای ‏ میله‌ای چنان‌چه مسیر سس یک 
مسیر خرابی باشد آن‌گاه مسیر زب‌ابزیا هر ترکیب 
دیگر این سه عضو نیز باعث خرابی سازه می‌شود. 
به‌عبارت دیگر عدم حضور سه عضو مذکور به‌طور 
همزمان (مستقل از ترتیب وقوع آن‌ها) منجر به 
ناپایداری خرپا خواهد شد. حال اگر در این خرپا دو 
عضو و 6 حذف شده باشد. قطعاً حذف عضو سوم 
() باعث خرابی است و نیازی به تحلیل خرپاو 
کنترل شرط خرابی نیست. همان گونه که بیان شد شرط 
خرابی خرپا در تحلیل به روش سختی. صفر شدن 
دترمینان ماتریس سختی کل و در تحلیل با روش نرمی 
پیشنهادی ایجاد نیرو در میله‌های حذف‌شده است. 

در اینجا ذکر دو نکته ضروری است. اول این کنه 
برای تعیین احتمال خرابی مسیر [ حس] حس1 باید 
خرپایی تحلیل شود که ابتدا عضو ام و سپس عضو > 
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عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


ام آن حذف شده است و نیازی به حذف عضو زام 
نیست. ثانیاً در یک سازه با 5 درجه نامعینی. هر مسیر 
شامل 5+1 عضو قطعاً یک مسیر خرابی است و نیازی 
به کنترل شرط خرابی يا استفاده از تکنیک‌های هوش 
مصنوعی نیست. 

در این تحقیق برای ذخیره‌سازی مسیر خرابی 
1 جح در خرپای 0 میله‌ای از رشته‌ای به طول 2 
استفاهه مت نود کته تسام کاراگزهتای: مشتاطر سا 
میله‌های یاد شده 0 و مابقی 1 است. 
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تعداد مسیرهای خرابی که با حذف اعضای با, 1 
به وجود می‌آید برابر با ۲۱-۲ است. ۵ مسیر دیگر که 
در ادامه‌ی عملیات شاخه‌ای به‌دست ی ایشند نیز با 
همین شکل تشکیل می‌شود و مورد شناسایی عامل 
هوشمند قرار می‌گیرد. بنابراین دیگر نیازی به تحلیل و 
ذخیره‌سازی تکراری این مسیرها نیست و در حجم و 
زمان محاسبات صرفه‌جویی می‌شود. به‌طور کلی هر 
رشته‌ی با طول # و دارای مد کاراکتر صفر, نشانگر 
یک خرپای 7 عضوی ناپایدار است که ,1 عضو خحود 
زاآن دشت فده انتت: اوق ترالن مخلفت: یات اسشاه 
رشته‌ی فوق می‌شود اما فقط توالی اول. تحلیل و 
ذخیره می‌شود و 0,۳1 توالی دیگر پس از جستجو در 
بانک خرپاهای ناپایدار شناسایی می‌شود و از تحلیل و 
ذخیره‌سازی تکراری آن‌ها خود داری می‌گردد. 

نحوه‌ی جستجو در بانک تأثیر عمده‌ای در زمان 
محاسبات دارد. در یک جستجوی ترتیبی باید اعضای 
پانک از اول به‌ترتیب با عضو جدید مقایسه شود تا در 
صورت یکسان بودن رشته‌ها جستجو پایان یابد. 
حداکثر تعداد عملیات لازم در این روش جستجو برابر 
با تعداد اعضای بانک است. تعداد عملیات مقایسه با 
جلو رفتن عملیات شاخه‌ای و زیاد شدن تعداد اعضای 
بانک افزایش می‌یاند. و اتتلاف قابل‌تنوجهی دز زسان 
به‌وجود می‌آید. به‌جای این روش می‌توان از روش 
سریع جستجوی دودویی استفاده کرد [32] لیکن باید 


سال بیست و چهارم. شماره دی ۱۳۹۲ 


قبلاً اعضای بانک به‌صورت صعودی یا نزولی مرتب 
شده باشد. با توجه به ذخیره‌سازی رشته‌ها به‌صورت 
بیان‌شده امکان مقایسه و مرتب‌سازی آن‌ها وجود دارد. 
اوق ها شا نمی شیاه دراه ۳ 
است. چنان‌چه بانک به‌صورت نزولی مرتب شده باشد 
برای پیدا کردن یک رشته باید آن را با عضو وسط 
بانک مقایسه کرد. اکنون اگر رشته‌ی مورد نظر از 
رشته‌ی وسطی بزرگ‌تر باشد جواب در نیمه‌ی بالایی 
بانک است و برعکس. با این روش در گام اول فضای 
جستجو به نصف آن و در گام دوم به یک چهارم تقلیل 
وتات ای رها تفای که تفه معط مد انیم 
ادامه می‌یابد. حداکثر تعداد عملیات لازم سم در این 
روش برای بانک با تعداد اعضای لا از رابطه‌ی )۱٩(‏ 
به‌دست می‌آید. 

۱۹( 1+[ رع10)]10ح م1 


در رابطه‌ی فوق 101 تابع جزء صحیح است. ور 
توا مه ابا مت انب کرک انس لا اناد 
از این روش در یسک فضای جستجو بسزرگ 
اجتناب‌ناپذیر است. مرتب‌سازی رشته‌ها در بانک از 
ابتدا باید به‌گونه‌ای باشد که هر رشته در جای خحود 
ذخیره شود تا نیازی به مرتب‌سازی بانک پس از ورود 
هر رشته جدید وجود نداشته باشد. در این تحقیق برای 
پیدا کردن موقعیت صحیح یک رشته‌ی جدید. از روش 
مرتب‌سازی فوری استفاده شده است [32]. 

همواره در حین انجام عملیات. زمان قابل توجهی 
ی کش زیخ کر نیز فی نوی با اساده از 
تکنیک‌های هوش مصنوعی می‌توان مسیرهای تکراری 
را شناسایی و از انجام عملیات تکراری خودداری 
نمود. این کار باعث می‌شود که مدت زمان مورد نیاز 
برای یافتن احتمال خرابی سازه کاهش چشم گیری یابد. 


عامل هو شمند فاز دوم 


در بهینه‌یابی با الگوریتم ژنتیک جمعیت اولیه به‌صورت 
کافلا تضادفی ترلید می‌ شود برای: تولید تسل‌هاق 


نشریه مهندسی عمران فردوسی 


بعدی رقابت بین کروموزوم‌ها صورت می‌پذیرد و 
کروموزوم‌های شایسته تر شانس بیش‌تری برای ادغام 
پیدا می کنند. به‌این‌ترتیب. در ادامه و پس از تولید 
نسل‌های بعدی کروموزوم‌های به‌تری ایجاد می‌شود. 
با دقت در این کروموزوم‌ها مشاهده می‌شود که شباهت 
بسیاری بین آن‌ها وجود دارد و حتی بسیاری از آن‌ها 
تکراری هستند. این کروموزوم‌ها يقیناً فرزنندان 
همان‌هایی هستند که با یب کون بالاتری داشته‌اند. 
در نسل‌های بعدی نیز دوباره تولید می‌شوند. 

اش مساله کر تسا هام: بایان سارز یر ااشت؛ 
به‌طوری که تعداد کپی‌های یک کروموزوم حتی از 
نصف تعداد جمعیت آن نسل بیش‌تر می‌شود. انجام 
عملیات تکراری برای یک کروموزوم خاص در یک 
نسل و نیز در نسل‌های بعدی هیچ فایده‌ای ندارد و 
فقط زمان برنامه را به هدر می‌دهد. 

دن انرت ین از غاما فوقسد برای یناب این 
کروموزوم‌ها استفاده شده است. عامل هو شمند فاز دوم 
در هر نسل کروموزوم‌هایی را که شایستگی آن‌ها از 
متوسط اش کی یبای فان مس نس اشیتت: ذخیره 
می کند. در صورت ترا سم کروموزوم‌ها عامل 
هو شمند آن‌ها را به‌ساد گی تشخیص می‌دهد و به‌جای 
تحلیل‌های اضافی و وقت‌گین مقدار شایستگی نظیر 
آن‌ها را مستقیماً اعلام می‌کند. 

منطق استفاده‌شده در این عامل هوشمند در عین 
سادگی, کاربرد بسیاری دارد. تنها نکته ایین است که 
برای ذخیره‌سازی رشته‌ها باید از الگوریتم‌های مناسب 
مانند آنچه در فاز اول معرفی شد. استفاده شود تا 
کارایی این عامل بالاتر برود. 

عامل هوشمند فاز دوم صرفا برای مسائل بهینه‌یابی 
تحت قید قابلیت افتماد کاربرد نذارد: از ایس عامل 
می‌توان براش تام سا ناس کته زمان تصلی از 
عملگرهای الگوریتم ژنتیک به‌طور قابل توجه بیش‌تر 
است استفاده کواگ: 


نشریه مهندسی عمران فردوسی 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


همان گونه که بیان شد در اغلب موارد مهندسی 
تحلیل و تعیین شایستگی کروموزوم‌ها حدود ۸۵ تا ۹۵ 
درصد زمان برنامه را به خود اختصاص می‌دهد. این 
عدد در مورد مسأله‌ی مورد بحث در اين تحقیق حتی 
از ۹٩‏ درصد بیش‌تر است؛ یعنی زمان لازم برای 
عملگرهای الگوریتم ژنتیک از یک درصد کم‌تر است. 
حال اگر از این یک درصد زمان صرف نظر شود 
نتیجه‌ی جالبی حاصل می‌گردد. در یک الگوریتم ژنتیک 
با تعداد جمعیت ۲و 6 نسل هر یک بار که عامل 
هوشمند کروموزومی تکراری پیدا می‌کند از محاسبات 
مجدد در مورد آن جلوگیری می‌کند. لذا به مقدار 
(1/)۳۵ در زمان محاسبات صرفه‌جویی می‌شود. 

در ادامهء مطابق مطالبی که در راهکار چهارم ارائه 
خواهد شد. مشخص می‌شود که مقدار صرفه‌جویی 
به‌وجود آمده بیش‌تر از این مقدار است. بر اساس نتایج 
بخش یادشده زمان تحلیل کروموزوم‌های با شایستگی 
بالا نسبت به سایرین به مراتب بالاتر است. به‌عبارت 
دیگر تحلیل کروموزوم‌هایی که شایستگی کم‌تری دارند 
خیلی سریح‌تر انجام می‌شود. 

مجددا یاداوری می‌شود که عامل هوشمند تنها 
شایسته‌ترین کروموزوم‌ها را ذخیره و شناسایی می‌کند؛ 
یعنی همان کروموزوم‌هایی که زمان تحلیل بالاتری 
دارند. بنابراین» صرفه‌جویی به‌وجود آمده در اثر 
شناسایی یک کروموزوم و عدم انجام محاسبات 
تکراری به‌مراتب بیش‌تر از مقدار (1/)۳0 است. 

مطالب عنوان‌شده با ذکر متال عددی شرح داده 
می‌شود. در یک مسأله‌ی بهینه‌یابی با جمعیت ۱۰۰ 
کروموزوم و ۱۰۰۰ نسل. چنان‌چه ۵۰۰۰ کروموزوم 
توسط عامل هوشمند شناخته و عملیات آن‌ها کنار 
گذاشته شود. تصور می‌شود که زمان محاسبات ۵۰ 
درصد کاهش يافته است. لیکن به دلیل آن‌که تحلیل 
این ۵۰۰۰ کروموزوم تکراری و شایسته بیش‌تر از 
سایرین طول می‌کشد زمان صرفه‌جویی شلده بیش‌تر 


خواهد بود. 
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عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


راهکار پیشنهادی در بهینه‌یابی 


در برخی از مسائل بهینه‌یابی در زمینه‌ی مهندسی در 
حین محاسبه شایستگی هر کروموزوم می‌شود و قبل از 
اتمام تحلیل می‌توان در مورد شایسته بودن آن 
کروموزوم قضاوت کرد. بسیاری از کروموزوم‌ها از 
اوایل تحلیل نشانه‌هایی از نقض قید از خود نشان 
می‌دهند. اگر مقادیر این نقض قیدها زیاد باشد, لزومی 
به اتمام تحلیل و محاسبه‌ی تابع پنالتی و غیره نیست. 
به‌عبارت دیگر در همین مرحله مشخص است که این 
کروموزوم شایستگی انتخاب برای ادغام را ندارد. 

از سوی دیگر برخی کروموزوم‌ها معرف طرح‌های 
دست بالا هستند به‌نحوی که تابع همدف را بهینه 
نمی کنند. این گونه کروموزوم‌ها معمولاً شایستگی 
پایینی دارند و شناسایی آن‌ها پیش از اتمام تحلیل در 
برخی مسائل مهندسی امکان‌پذیر است. اتمام تحلیل 
این گونه کروموزوم‌ها عمللا فایده‌ای ندارد و فقط زمان 
محاسبات را بالا می‌برد. 

هسدف اصلی راهکار چهارم اعمال چنین 
قضاوت‌هایی است که نتبجه‌ی آن جلوگیری از 
محاسبات اضافی و صرفه‌جویی در زمان خواهد بود. 
خوشبختانه روابط تعیین احتمال خرایبی یک سیستم 
سازه‌ای به‌نحوی است که می‌توان در میانه (يا حتی در 
اوایل) تحلیل در مورد وضعیت سازه اظهار نظر نمود. 

و 
قرار گرفت. مطابق مطالب بیان‌شده در آن بخش شرط 
لازم برای تولید مسیرها (عملیات شاخه‌ای) این است 
که احتمال وقوع مسیر مورد نظر از ۷ ۸/۵۷۳,10 
کوچک‌تر نباشد. ,18*۳ حداکثر احتمال مسیرهای 
خرابی گذشته است. 

در هر مرحله از عملیات شاخه و کرانه با یافتن 
مسیری محتمل تر» مقدار احتمال آن در ,12۳ قرار 
داده می‌شود. به‌این‌ترتیب مقدار ۲ ۷2۳,10 نیز 
افزایش می‌یابد. در این‌جا باید دقت شود که برخی از 
شیرهایی کر اشدای عملانت: شتا هی کرانه یشان 


شاخته شده و تولید شده‌اند. از مقدار حدید 


سال بیست و چهارم. شماره دو. ۱۳۹۲ 


2۳,7 کوچک تراند و طبق عملیات کرانه‌ای 
نبایستی تولید می‌شدند. 

از این مطالب چنین استنباط می‌شود که اگر مقدار 
۷۵7 که در پایان عملیات به‌دست می‌آید ( ۳ از 
ابتدا معلوم بود. تعداد مسیرهای حذف‌شله بیشتر 
می‌شد. یعنی بسیاری از مسیرها بدون این‌که تأثیر 
یو سل دقن اسقیارد اتسیینال ام اش دافته ماش فرلید 
شده‌اند و در نتيجه زمان زیادی به هدر رفته است. 

قو بهیصیاییساژه‌ها با انففاده او الکو ریت وتتینک 
جمعیت اولیه ب‌صورت کاملاً تصادفی تولید می‌شود. 
برخی از این کروموزوم‌ها به‌شدت قید قابلیت اعتماد را 
نقض می‌کنند. در حالی که برخی دیگر احتمال خرابی 
بسیار کوچکی دارند. 

راهکار چهارم پیشنهادی با یک مثال شرح داده 
می‌شود. یک خرپای ۸ میله‌ای نامعین مفروض است. 
همدف یافتن سطح مقطع میله‌های خرپااست 
به‌ طوری که کم‌ترین وزن برای خرپا حاصل شود و 
مقدار احتمال خرابی کل آن از * ۱۰ (احتمال خرابی 
مجاز) کوچک‌تر باشد. مقدار 1 برابر با ۳ در نظر 
گرفته می‌شود. 

طبق مطالب بیان شده در هر یک از کروموزوم‌ها 
در طول روند بهینه یابی. تمام مسیرهایی که احتمال 
اوق ای ۱۳ کتکایر استاناید کار گذاشقه کتود 
اکنون چنان‌چه در یکی از کروموزوم‌ها(خرپاها) 
احتمال خرابی محتمل‌ترین عضو کوچک‌تر از ۱۰۲ 
باشد (مثلاً ۲ ۱۰۳) احتمال وقوع تمامی مسیرهای 
خرابی قطعاً از عدد مذکور کوچک‌تر است. طبق 
عملیات کرانه‌ای کلیه‌ی این مسیرها در خرپای مورد 
نظر پیش از تولید باید کنار گذاشته شوند. عدم تولید 
این مسیرها زمان محاسبات را به‌طور قابل‌توجهی 
کاهش می‌دهد. 

به‌طور خلاصه معیار حذف مسیرها در راهکار 
چهارم برآورده کردن رابطه‌ی (۲۰) است. 


۲ ۱۷۵/0 ,, ,0۵:6۳ ک و۲ 


211۰ 


۳ 
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در رابطه‌ی فوق ,,:۲ احتمال خرابی مجاز (قید 
مسأله‌ی طراحی) است. ,۷1۵ و بر در عملیات شاخه 
و کرانه شرح داده شد. این معیار برای کروموزوم‌هایی 
که سطح مقطع میله های آن‌ها دست بالا است کاربرد 
داره. بتلایهی استنا کسه.میران جریمه‌ی این گونه 
کروموزوم‌ها صفر است لیکن شایستگی مناسبی ندارند. 

هنوز امکان صرفه‌جویی بیش‌تر در حجم و زمان 
محاسبات وجود دارد. در مشال بادشده تمامی 
کروموزوم‌هایی (خرپاهایی) که احتمال خرابی آن 
بزرگاتر از ۱:7 اسنت نقض قید کنرده و یایند جریمه 
شوند. لیکن بین دو کروموزوم که احتمال خرابی آن‌ها 
۱/3 و ۰/۹۹ است تفاوت زیادی وجود دارد. 
هم‌چنین میزان نقض قید و جریمه‌ی این دو کروموزوم 
کاملاً متفاوت است. در هنگام تحلیل کروموزوم اول تا 
لحظه‌ی پایان نمی‌توان پیش‌بینی نمود که آیا نقض قید 
می‌کند یانه. اما در مورد کروموزوم دوم وضعیت 
متفاوت است. اگر مجموع احتمال مسیرهای خرابی 
(کران فوقانی کرنل) در اوایل تحلیل به عددی مانند 
برسد. نیازی به ادامه‌ی تحلیل نیست و در همین 
مرحله مشخص است که نقض قید صورت پذیرفته و 
چنین کروموزومی باید جریمه شود. 

این امر سبب می‌شود که در بسیاری از موارد 
عملیات زودتر از موعد مقرر پایان یابد و بخش عظیمی 
از مسیرهای خرابی و محاسبات نظیر آن‌ها کنار گذاشته 
شود. تنها نکته‌ی باقیمانده‌ی تعریف معیاری دقیق برای 
اتمام عملیات است. از رابطه‌ی (۲۱) می‌توان برای 
کنترل شرط پایان استفاده نمود. 


)۳۱( طبار 
1<1 


عبارت سمت چپ نامعادله نشان دهنده‌ی مجموع 
احتمالات مسیرهای خرابی حاکم از ابتدای تحلیل تا 
مسیر جاری ۳ است. ,۳ احتمال خرابی مجاز (قید 
مسأله‌ی طراحی) و یک ثابت عددی است که 
انتخاب آن به‌عهده‌ی طراح است. کروموزومی که در 
حین تحلیل معیار رابطه‌ی (۲۱) را ارضا می‌کند باید 
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بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


به‌شدت جریمه شود به‌طوری که شانس انتخاب و ادغام 
پیدا نکند. انتخاب عددی مانند ۲ با ۳ برای » انتخاب 
تام اش زوا کروشو ای که با ایب ان 
رابطه (۲۱) را ارضا کند. یقیناً کروموزوم شایسته‌ای 
نیست. معیار رابطه‌ی فوق‌الذکر مختص کروموزوم‌هایی 
است که سطح مقطع میله‌های آن‌ها بسیار ضعیف است. 

به‌طور خلاصه با استفاده از راهکار چهارم. هر چه 
احتمال خرابی یک کروموزوم از مقدار احتمال خرابی 
مجاز (قید مسأله) دورتر باشد به تحلیل کم‌تری نیاز 
دارد. بنابراین راهکار اول در نسل‌های اول الگوریتم 
ژنتیک اثر بیش‌تری دارد. 


روند حل مسأله و ارتباط بین راهکارهای 
پیشنهادی 

در این قسمت روند کلی حل مسأله بهینه‌یابی و ارتباط 
بین راهکارهای پیشنهادی به‌طور خلاصه به‌شرح زیر 
ارائه می‌شود. 

در بهینه‌یابی به‌وسیله‌ی الگوریتم ژنتیک تولید یک 
نسل جدید منوط به تعیین شایستگی کروموزوم‌ها (در 
این تحقیق خرپاها) و سپس انتخاب شایسته‌ترین‌ها 
انست::دن مساله‌ی؛ حاضی شانستگی غریاها بیستگی یه 
مقدار احتمال خرابی سیستم سازه‌ای دارد. همان گونه 
که بیان شد تعیین احتمال مذکور به‌طور دقیق تقریباً 
مستلزم صرف زمان بسیار زیاد و تقریبا غیر عملی 
است. بنابراین در تحقیقات گذشته از فرضیات بسیار 
ساده‌کننده و محافظه‌کارانه برای تعیین احتمال خرابی 
استفاده ده استخ: ,ورين یی با اشفاده از زاهک ار 
اول تعداد کثیری از مسیرهای همبسته حذف می‌شوند 
و احتمال خرابی سیستم به‌طور سریع‌تر از روی احتمال 
وقرع مسیرهای خرابی حاکم تعیین می گردد. برای 
تشکیل مسیرهای خرابی باید خرپاها را به‌کرات تحلیل 
سازه نمود. 

راهکار دوم اين تحقیق امر تحلیل سازه‌ی خرپا را 
تسریع می‌کند. استفاده از این راهکار تأثیری بر مقدار 


سال بیست و چهارم شماره دو» ۱۳۹۲ 


عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


احتمال خرابی ندارد و تنها زمان محاسبات را به‌طور 
قابل توجهی کاهش می‌دهد. 

راهکار سوم نیز تأثیری در جواب‌های به‌دست 
آمده ندارد و تتهااز تشکیل مسیرهای خرابی و 
کروموزوم‌های تکراری و محاسبات اضافی نظیر آن‌ها 
جلوگیری می‌کند. 

با استفاده از راهمکار چهارم تعداد زیادی از 
گرزموزوم‌هایی که شایستگی مناسبی تدازشنه در تحین 
لیا ۵ هنارای قطان خوانیای سیر تیان 
محاسبات بیهوده جلوگیری می‌گردد. 


مثال‌های عددی 


در همه‌ی مثال‌های ارائه‌شده به‌جای احتمال خرابی 
تتاژه او کزان فرقای کریل ابرفاده شوه استه 

پیش از ارائه‌ی مثال جامع بهینه یابی ابتدا یک 
مثال برای ارزیابی راهکارهای اول تا سوم مطرح 
می‌شود. در این مثال با استفاده از راهکارهای بادشده 
احتمال خرابی یک خرپای ۲۵ میله‌ای تعیین می‌شود و 
سهم هر یک از راهکارها مشخص می‌گردد. پس از آن 
در مثال دوم تأثیر راهکار چهارم در خرپای ۲۵ میله‌ای 
مثال اول بررسی می‌شود. سپس با توجه به نتایج 
حاصل و بااستفاده از راهکارهای پیشنهادی مثال بهینه 
یابی خرپای یادشده تحت قبد قابلیت اعتماد بررسی 
می‌شود. در مثال بهینه‌یابی متغیرهای طراحی سطح 
مقطع میله های خرپا و قید مسأله کوچک‌تر بودن 
احتمال خرابی سازه از مقدار احتمال خرابی مجاز 
است. در همه‌ی مثال‌ها (مطابق مطالب یادشده در 
بخش دوم) متغیرهای تصادفی عبارتند از بارهای وارد 
بر سازه و مقاومت میله های خرپا. کلیه متغیرهای 
تصادفی نرمال و غیرهمبسته فرض می‌شوند. 


ال ۱. خرپای ۲۵ میله‌ای مطابق شکل (4) مفروضص 
است. سطح مقطع اعضاء مختصات گره‌ها و میانگین 


بارهای وارد به‌ترئتیب در جداول (۲-۵) ارائه شده 


سال بیست و چهارم شماره دو, ۱۳۹۲ 


امستام یب بر آکنند کین بانها ۰۲۲ انشتهه میا کین و 
ضریب پراکندگی تنش تسلیم اعضا به ترتیب برابر با 
۰ کیلوگرم بر سانتیمتر مربع و ۰/۰۳ است. تنش 
تسلیم اعضا و بارهای وارد متغیرهای تصادفی نرمال 
غیرهمبسته و رفتار مصالح در کشش و فشار یکسان 


شکل ۶ خرپای ۲۵ میله‌ای مثال ۳ 


جدول ۳ سطح مقطع اعضای خرپای مثال ۳ 


شماره‌ی عضو سطح مقطع (سانتی‌متر مربع) 
۱ ۳/۸۹ 
۳4۵ ۳/۳/۳۳ 
۳ ۲ 
14 ۳/۳۵ 
۷4۸ 9۸۱۷۹ 
۳ ۱۳۸۵ 
۱۳۵۳ ۱/۰ 
۱:۵۷ ۳/۳ 
او( ۱/1۵ 
۱۱۸۵۱ ۰/۸۶ 
۱۹۲۰ ۱/۳۸ 
۳۳۵ ۱۳/۰۲ 
۳۳۵۶ 1/34 


مقدار همبستگی شاخص فشردگی با مشخصه‌های 
مرتبط با ساختمان اولیه‌ی نمونه‌ها به‌مراتب بیش‌تر از 
مقدار همبستگی آن با مشخصات مرتبط با جنس ذرات 
نظیر خواص خمیری خاک است. به عبارت دیگر 
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روابط مبتنی بر شرایط ساختمان اولیه نظیر نسبت 
پوکی. دانسیته‌ی خشک و رطوبت که همگی به‌نوعی 
بیانگر فاصله‌ی ذرات خاک از یکدیگر و نهایتاً پتانسیل 
نشست هستند در تخمین مقدار شاخحص فشردگی 
کارایی بیش‌تری دارند. هم‌چنین روابطی که تنها مبتنی 
بر جنس و خواص ذرات ریزدانه (نظیر حدود اتربرگ) 
و بدون توجه به تخلخل اولیه‌ی آن‌ها می‌باشند. از نطر 
آماری معنی‌دار نیستند و قادر به تخمین درست 


شاحص فاد کی نمی‌باشند. 


حدول ۶ مختصات گره‌ها در مثال ۱ سانتی‌متر) 


شماره‌ی گره ۹ ۷ 7 
۱ ۰ 9/۲۵ 9۸ 
۲ ۰ ۵/- 9۸ 
۳ ۹۰9/۳۵ ۹۵/۳۵ ۳0 
3 ۹9/۹6 9/۵- ۳۶ 
۵ ۵/- ۵/- ۳۹۶ 
. ۵/- 9/۲۵ ۳۹ 

۰ ۳۶ ۳ ۷ 

۰ "۳ ۳۶ ۸ 

۹ ع- 5" ۰ 

۰ ۳۶ "8 ۲۰ 


جدول ۵ داده های آماری بارهای وارده مثال ۱ (کیلوگرم) 


شماره‌ی گره | جهت >« جهت ۲ جهت 7 
۱ ۸۳۹۰ ۸۳۹۰ ۳۳۹۰" 
۲ ۰ ۰ ۳۳ 
۳ ۸۱۸۹۹۰ ۸۱۸۹۹۰ 
۵ ۰ ۰ 


تحلیل این سازه بدون استفاده از راهکار ۱ حدود 


۳ ساعت طول می‌کشد و با جمع چندین میلیون مسیر 
خرابی که اکثبر آن‌ها همبسته اسست» عددی بسیار 


محافظه‌کارانه برای کران بالای احتمال خرابی سیستم 
سازه‌ای نعیین می شود. در صورت استفاده از راهکار ۱ 
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بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


(با کاهش حدود ۰ درصد در زمان محاسبات از ۳ 


ات لمات 


جدول 1 تأثیر راهکارهای پیشنهادی در مثال ۱ 


راهکار آراهکار ۱|راهکار ۱] با هر سه 
۱ ۳ و۲ | راهکار 
زمان محاسبات (ثانیه) | ۲۱/۱ | ۱۲/۷۰ | ۲۰/۰۱ 1 ۱۱/۵۶ 
محتمل 
تعداد تحلیل قابلیت 


۲۰ ۲۰ ۲۰ ۳۰ 


۶ | ۱۶۹۶ ۱۸۶ ۱۱۸۶ 
اعتماد 


با ملاحظه‌ی نتایج به‌دست آمده از جدول (0) 
مشاهده می‌شود که استفاده از راهکارهای ۰۲ ۳و 
ترکیب آن‌ها باعث کاهش چشم‌گیر زمان محاسبات 
سل ۵ انلنتتا: بیش‌ترین مان وت ی رت فادها 
راهکار ۱ به‌دست می‌آید. در صورت استفاده از راهکار 
۲ یعنی بهبود روش جبری نیروهاء بدون این‌که تغییری 
در تعداد تحلیل‌ها ایجاد شود زمان محاسبات کاهش 
یافته است. صرفه‌جویی به‌وجودآمده در اثر استفاده از 
عامل هوشمند و جلوگیری از محاسبات اضافی است. 
مراجع [5 ,2] احتمال خرابی سیستم سازه‌ای این خرپا 
با ۱/۰۰۱۰ تعیین کرده‌اند. کران فوقانی احتمال 
۳ طای چ 4 2 
خرابی در این تحقیق برابر با ۰ تعیین شد که 
نسبت به تخمین محافظه کارانه‌ی مراجع یادشده بهیود 
یافته است. 
مثال ۳. هدف از ارائه‌ی این مثال بررسی تأثیر راهکار 
چهارم یعنی استفاده از قید احتمال خرابی مجاز در 
زمان محاسبات است. بدین منظور خرپای ۲۵ میله‌ای 


مثال (۱) با سه دسته سطح مقطع متفاوت مورد تحلیل 
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عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


قرار می گیرد. سطح مقطع خرپای (الف) در جدول (۳) 
قید شده است. خرپای (ب) دارای میله‌های قوی با 
سطح مقطع ۱۰ سانتی‌متر مربع و خرپای (ج) دارای 
میله‌های ضعیف با سطح مقطع ۳ سانتی‌متر مربع است. 
برای صرفه‌جویی در زمان از همه‌ی راهکارهای ۲۰۱ 
و ۳ (فاز اول) در تعیین احتمال خرابی استفاده می‌شود. 
مقدار 11 و !| به‌ترتیب برابر با ۲ و ۳ در نظر گرفته 


می‌شود. 
جدول ۷ بررسی تأثیر راهکار ؛ 
بدون راهکار ۶ | با راهکار ۶ 
زمان محاسبات ۱۱/۶ ۱۰/۸۹ 
خرپای الف 
۱ تعداد مسیرها ۳۰ 1۹0 
زمان محاسبات ۰/۱ ۱/۰ 
خرپای ب 
تعداد مسبر‌ها ۱۰ ۰ 
زمان محاسبات | ۳ ۷/۸۱/۱۰ ۲ ۶/۷۷/۱۰ 
ری کج 
تعداد مسیرها 3 ۱ 


با مشاهده‌ی جدول (۷) اولین نتیجه‌ای که حاصل 
می‌شود این است که استفاده از راهکار چهارم زمان 
محاسبات را کاهش می‌دهد. نتیجه‌ی دوم این است که 
زمان تخطیل دز عرپا با سطم نقطم :تفاوت با یکنذیگر 
فرق دارد زیرا مسیرهای خرابی و تعداد مسیرها در دو 
خرپا با یکدیگر متفاوت است. با اندکی تغییر در سطح 
مق یلها تیان طمیرمای مان ماش بات بط و 
چشم‌گیر تغییر می‌کند. این امر در بهینه‌یابی خرپاها 
تحت قید احتمال خرابی تأثیر قابل توجهی دارد. 
نتیجه‌ی سوم اين که تأثیر راهکار چهارم روی دو 
خرپا با سطح مقطع متفاوت با یکدیگر فرق دارد. دلیل 
این امر در بخش 1-4 بیان شد. در مورد خرپای (ب) 
که دست بالا طراحی شده است. پس از اولین تحلیل 
مشخص می‌شود که احتمال خرابی محتمل‌ترین عضو 
۳ است. تمامی مسیرهای خرابی که در 
ادامه‌ی عملیات شاخه‌ای حاصل می‌شوند احتمال 
کوچک تری دارند لذا توسط راهکار چهارم از تولید 


آن‌ها جلوگیری می‌شود. 


سال بیست و چهارم شماره دی ۱۳۹۲ 


در مورد خرپای (ج) پس از تشکیل اولین مسیر 
خرابی. احتمال آن برابر با ۰/۳۷ تعیین می‌شود. با توجه 
به مفاهیم راهکار چهارم دیگر نیازی به ادامه‌ی تحلیل 
نیست. زیرا در همین مرحله ناشایسته بودن این خرپا 
معلوم می‌شود. 

در مورد خرپای (الف) با استفاده از مفاهیم راهکار 
چهارم مسیرهای خرابی بیش‌تری کنار گذاشته شده 
است. 
مشال ۳ :در این متال با استفاده از راهکارهنای 
پیشنهادی به بهینه‌یابی خرپای ۲۵ میله‌ای مثال ۱ تحت 
قید قابلیت اعتماد سیستم سازه‌ای پرداخته می‌شود. 
احتمال خرابی مجاز (قید مسأله) برابر با ۱/۱۰ 
است. نتایج بسب‌دست: آمته در جدول (۸) با مراجع 
[5 ,2] مقایسه شده است. در مراجع یادشده احتمال 
خرابی سیستم سازه‌ای در جهت اطمینان برابر با 
مجموع احتمال خرابی اعضا فرض شده است. 


جدول ۸ نتایج بهینه‌یابی مثال ۳ و مقایسه با مراجع دیگر 


سطح مقطع (سانتی‌متر مربع) 
۱ 7۷-۷۲ , 0.05 -0.20 
مرجع [2]۳ مرجع [۵] 0 مرجع [۵] این تحقیق 

32 1۱۳۸۷ 3/۹۰ ۳۹ ۱ 
۳۰ 3 ۸۵ 3 ۲۵ 
2/۰ ۱۷/۵۰ 3/۹۰ ۱۷/۷ ۳۶ 
۹۰ 1۱۳۳۹ 3/4 ۳/۳۹ 14 
۹۰ ۱/۹۹ ۹/۹9 ۱/0۲ ۷/4۸ 
۳/۸ ۱/۳ ۹۰ ۱۳۸ ۱۱۲۰۰۱ 
۳/۵۸ ۱/۷۸ ۳/۸ ۱۷۳۷۷ ۱۳۳ 
۳/۳۰ 1۱۰۰ ۳/۸ 1/۸۸ ۱:0۷ 
۳/۳۰ ۳/۱۱۷۹ 1/۳۱ ۱/۳/۳۹ ۱۳۹ 
۱/۳/۳۹ ۱۱/۳/۳۰ ۳/۷۹ ۱/۷۸ ۱۱۸۰۲۱ 
۱/۳/۳۹ ۳/۵0۹۵ ۳/۳۹ ۳/۳۹۳ ۱۹۳۰ 
1۳۸ 1۹۳ 3/۹۰ #۹ ۳۳۵6 
1۳۸ ۱۷/۸۳ ۹/۹ ۷/۳۹۹ ۳۳۵ 
٩۱/۹۵ | ۸۱ ۱۱۸/۷ ۱ ٩۷۸ | )162( وزن‎ 


نشریه مهندسی عمران فردوسی 


با توجه به نتایج جدول (۸) مشاهده می‌شود که در 
این تحقیق طرح بهینه‌ی مطلوب‌تری نسبت به مراجع 
[5 ,2] به‌دست امده است. دلیل این امر تخمین به‌تر 
نسبت به این مراجع برای احتمال خرابی سیستم 
سازه‌ای و استفاده از الگوریتم ژنتیک جزیره‌ای است. 

شایان ذکر اسبت کته معفرهسای طراضی در این 
تحقیق از نوع گسسته و در مراجع [5 ,2] از نوع 
گسسته و پیوسته است. جواب‌های ب‌دست آمده در 
هر صورت از هر دو حالت مراجع یادشده به‌تر است. 

اکنون برای ارزیابی راهکارهای پیشنهادی و سهم 
آن‌ها در تسریع عملیات بهینه‌یابی زمان محاسبات ۵۰۰ 
نسل الگوریتم ژنتیک با استفاده از هر یک در جدول 
(۵) ارانه می شون 

ذکر دو نکته در مورد اعداد این جدول ضروری 
است. نکته‌ی اول اين‌که زمان‌های درج‌شده در ایین 
جدول مربوط به کندترین پردازشگر (جزیره) است. 
علی‌رغم استفاده از پردازشگرهای یکسان (مطابق 
جدول ۱ با توجه به تفاوت کروموزوم‌های این 
جزایر زمان بهینه‌یابی در این جزایر متفاوت است. 

نکته دوم این که در دو بار بهینه‌یابی تحت قید 
قابلیت اعتماد سیستم سازه‌ای با شرایط کاملاً یکسان, 
ممکن است زمان محاسبات متفاوت باشد. دلیل این امر 
مجدداً تفاوت کروموزوم‌ها است. بنابراین زمان‌های 
قیدشده در جدول )٩(‏ قطعی نمی‌باشد. بلکه تنها 
معیاری برای مقایسه است. 


جدول ٩‏ تأثیر راهکارهای پیشنهادی در مثال ۳ 


راهکار زمان محاسبات (ثانیه) 
3 29۹9۵ 
:۶ ۳۳/۸ 
۱ الف. ۶ ۸۱۳/۳۲ 
۱ الف. ۶ 2 
۱ اف آب. ۶ ۷۸۷۸ 


تشربه مهنلسی خمران فردوسی 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


با توجه به نتایج جدول فوق استفاده از راهکارهای 
پیشنهادی باعت کاهش زمان محاسبات شده است. 
هرچند که سهم آن‌ها متفاوت است. 

در هن جز یره تعداد کل کروموزوم‌ها دی هر تسل 
برابر با ۲۰ در نظر گرفته شده‌است لذا تعداد تحلیل‌ها 
به‌ازای ۵۰۰ نسل برابر با :۱۰۰۰ می‌شود. این در 
حالی‌است که تعداد کروموزوم‌هایی که شایستگی آن‌ها 
توسط عامل هوشمند فاز دوم مورد ارزیابی قرار 
و باظری عیبر با ۵00۷ آیشت یی 
شایستگی بیش نیمی از کروموزوم‌ها توسط عامل 
هوشمند فاز سوم تعیین می‌شود و زمان محاسبات به 
نصف تقلیل می‌یابد. از طرفی میزان کاهش زمان 
وا من سم لاتم دار شاه یا کته 
کروم وزوم‌هایی توسط عامل هوشمند شناسایی 
می‌شوند. کروموزوم‌های شایسته هستند. با توجه به 
مطالب بیان‌شده در بخش 1-6 و نتیجه‌ی مثال (۲)» این 
کروموزوم‌ها زمان بسیار زیادتری نسبت به سایرین 
صرف می‌کنند. لذا عدم تحلیل این گونه کروموزوم‌ها 
زمان محاسبات را به حدود یک‌سوم کاهش می‌دهد. 
در شکل (۵) نحوه‌ی همکاری جزایر در یافتن جواب 
بهینه مربوط به جدول (۸) نمایش داده شده است. 
همان‌گونه که در شکل مشاهده می‌شود در مراحلی که 
انتقال بین جزایر صورت می‌پذیرد» جهش‌هایی در 
نمودارها به‌وجود آمده است. جواب بهینه در نسل ۲۸۸ 
در یکی از جزایر به‌دست آسده و در نسل‌های ۳۰۰ 
۰ و ۶۰۰ به جزایر دیگر انتقال داده شده است. 


1 51 101 181 201 23 301 3581 401 451 1 


۱۱۷۱۱۵۵۲ 0 ۹ 


شکل ۵ نحوه‌ی کارکرد جزایر موازی با یکدیگر 


سال بیست و چهارم شماره دو» ۱۳۹۲ 


عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


برای بررسی و مقایسه‌ی بیش‌تر نتایج بهینه‌یابی 
خرپای ۲۵ میله‌ای به‌ازای ضرایب پراکندگی مختلف 


بار و مقاومت در حدول (۱۰) ارائه شده استتا. 


جدول ۱۰ اثر ضرایب پراکندگی در خرپای منال ۳ 


شلع ویو رن ات 
عضو | مرجع [۵] | این تحقیق | مرجع [۵] | این تحقیق 
۱ 2/۱۸ ۷/۵ ۷۳/۹۸۳ ۸/۲ 

۸/۸ ۸۱/۸۳۳۶ ۷/۵ ۱/۷۹ ۳۵ 
۸/۸ ۱/۳۰۱ ۷/۵ ۹/۸۳۹ ۳ 
۱۰/۹ ۹/۳۶ ۱۰/۵ ۱۷/۸/۳۶ 14 
۱/۹ ۱۳/۹۶ ۱/۵ ۱/۳۳۵ ۷4۸ 
17 1۳۷۹ 1/۵ ۳۸۷۹۷ ۱ 
1/۰ ۲۳/۳۹ ۳/۹ ۲/۷۹ ۳ 
1/3 ۷۳/۸۷۲ 0/۷ 3/۲ ۷ 
1/۹ 11:۹ 0/۷ ۳/۳۱۸ ۰ 
۳/۷ 0۱۳۹ ۳/۹ ۳/۵۰۱ 
۳/۷ 1۳۵ ۳۹ 0/۱۹۸ ۰ 
۱/۲ ۱/۳۳۸ ۷/۵ ۹/6۰ ۱ ۵ 
۱/۲ ۱۳/۹/۷ ۷/۵ ۱۲/۰۱۹ ۳ ۶ 
۱9۸/۳ ۱۸۳/۸۹ ۱۳/۸۷۸ ۱۵۲۱۲ | 
)12( 


با توجه به جدول (۱۰) ملاحظه می‌شود که وزن 
خرپای بهینه در این تحقیق در تمامی حالات کاهش 
يافته است. هم‌چنین با نگاه کلی به نتایج این جداول و 
جدول (۸) مشاهده می‌شود که هر چه ضرایب 
پراکندگی بار و مقاومت بیش‌تر شوند. وزن خرپای 
بهینه افزایش می‌پابد. 
مشال. . در این مثال ابتدا تأثیر واهکارهای پیشنهادی 
در تحلیل قابلیت اعتماد گنبد ۱۲۰ میله‌ای (مطابق شکل 
۲ ارزیابی می‌شود و پس از آن به بهینه‌یابی این سازه 
با استفاده راهکارهای ارائه‌شده پرداخته می‌شود. ابعاد 
نشان داده‌شده در شکل مذکور بر حسب سانتی‌متر 
است. اعضای سازه مطابق شکل مذکور گروه‌بندی شده 
و سطح مقطع نظیر آن‌ها قن ول ار ۱۱ )واه کته 
است. تنش تسلیم اعضا و بارهای وارد متغیرهای 


تصادفی نرمال غیرهمبسته است و رفتار مصالح در 


سال بیست و چهارم» شماره دی ۱۳۹۲ 


۳۱ 


کشش و فشار یکسان فرض می‌شود. به کلیه‌ی نقاط 
آزاد این گنبد (نقاط غیر تکیه‌گاهی) بار قائم وارد 
می‌شود. میانگین بار قائم وارد بر تاج این گنبد برابر 7۰ 
کیلو نیوتن بارهای وارد به گره‌های تراز میانی ۳۰ کیلو 
نیوتن و بارهای وارد به گره‌های تراز تحتانی ۱۰ کیلو 
نیوتن فرض می‌شود. ضریب پراکندگی کلیه‌ی بارها 
است. میانگین و ضریب پراکندگی مقاومت تسلیم 
اعضا به‌ترتیب برابر با ۲۷۹۰ کیلوگرم بر سانتی‌متر مربع 
و ۰/۰۳ است. 

تحلیل این سازه نیز همانند خرپای ۲۵ میله‌ای 
دون اسفاده از راهمکار اول ساعت‌ها بسهطول 
می‌انجامد و از سوی دیگر جمع احتمال وقوع چندین 
میلیون مسیر خرابی که اکثر آن‌ها همبسته است» عددی 
بسیار محافظه‌کارانه برای کران بالای احتمال خرابی 
سیستم سازه‌ای نتیجه می‌دهد. در حالی که در صورت 
اسقاقه از واهکار ۱ مان میات یه ۲۸۲۷۸۹ کائتة 
کاهش می‌یابد. هم‌چنین کران فوقانی احتمال خرابی با 
استفاده از تکنیک‌های پیشسنهادی برابر با 
۰ محاسبه می‌شود که مرائب کوچک‌تر و 
به‌تر از مجموع احتمال خرابی اعضا یعنی ۳/۲۹۱۰ 
است. نتایج تحلیل قابلیت اعتماد این خرپا با استفاده از 
راهکارهای پیشنهادی دوم سوم در جدول ۱۲ ارائه 
شده است. اعداد این جدول بر اساس ۲ < 7 به‌دست 


ام اننتتتا. 


شکل ٩‏ کید ۱۲۰ میله‌ای ما ۶ 


نشریه مهندسی عمران فردوسی 


۳۲ 


حدول ۱۱ سطح مقطع اعضای گروه گنبد مثال ۶ 


شماره | سطح مقطع 
گروه | (سانتیمتر مربع) 
۱ ۳/۲۰ 
۲ 1۰ 
۳ ۳/۲۳۰ 
1/۸۰ 
۵ 1/۰۰ 
3 1/1۰ 
۷ ۳/۰۰ 


جدول ۱۲ تأثیر برحی راهکارهای پیشنهادی در مثال ۶ 


راهکار ۱ |راهکار ۱| راهکار ۱ | با هر سه 
و ۲ و ۳ راهکار 

زمان محاسبات 
۹ ( ۱۵۳ ۲۲۹۱/۸ ۱۱3۹/۰ 

(ثانیه) 

تعداد مسیرهای 

1۸۰ ۶۸۰ ۶۸۰ 1۸۰ 
تعداد تحلیا 

۱۳:۷ ۲۳۶۷ | ۷۳ ۳ ِ 

قابلیت اعتماد 


بارزترین نکته‌ای که از جدول (۱۲) استنباط 
شوه از اتف که اسیفادی ان راهکار مسیروت 
بهبود یافته‌ی نیروها تأثیر بسیاری بر زمان تحلیل داشته 
است. سازه‌ی گنبدی مورد بحث ۱۱۱ درجه آزادی 
فعال گرهی و تنها ٩‏ درجه نامعینی استاتیکی دارد. اینن 
بدان معناست که در روش سختی دستگاه معادلات 
ما دوه رصان که اه تون ری 
پیشنهادی تنها به حل معادلات 4 احتیاج است. این 
تفاوت فاحش با توجه به کثرت دفعات تحلیل تأثیر 
بسیار مطلوبی داشته و حدود ۹۵ درصد زمان برنامه 
(بدون هیچ تغییری در تعداد مسیرهای خرابی با در 
ت زه فانک 

تعداد تحلیل‌های قابلیت اعتماد با توجه عملکرد 
عامل هوشمند فاز اول از ۲۶۷۳ به ۲۳۶۷ کاهش یافته 


اشتتا: 


نشریه مهندسی عمران فردوسی 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


حال با استفاده از این راهکارها به بهینه‌یابی گنبد 
۰ میله‌ای تحت قید قابلیت اعتماد پرداخته می‌شود. 
احتمال خرابی مجاز (قید مسأله) برابر با ۱/۱۰ 
است. در مسأله‌ی بهینه‌یابی مقاومت میله‌ی متغیرهای 
تصادفی با توزیع لوگ نرمال و ضریب همبستگی برابر 
با ۰/۵ فرض می‌شود. هم‌چنین بارهای وارده به سازه 
همگی متغیرهای تصادفی کاملاً همبسته و با توزیع 
نرمال در نظر گرفته می‌شوند. 

هلف از ارات ای ال آفبات:عسی سرد آسر 
بهینه‌یابی خرپاهای بزرگ تحت قید قابلیت اعتماد 
سازه‌ای است که با استفاده از راهکارهای پیشنهادی 
ممکن می‌گردد. تعداد تحقیقات انجام‌شده (در زمینه‌ی 
بهینه‌یابی تحت قید احتمال خرابی سیستم سازه‌ای) که 
در آن‌ها ب‌صورت موفقیت‌آمیز به حل مثال‌های عددی 
بزرگ پرداخته شود کم بوده و تنها قاب‌ها و خرپاهای 
بسیار کوچک مورد بهینه‌سازی قرار گرفته است [33]. 
در اين مثال هم‌چنین سعی شده است از فرضیاتی 
واقعی‌تر در مورد نوع توزیع و همبستگی متغیرهای 
تصادفی استفاده شود. 

در این مثال از پنج بزدازافت‌کن موازی در امر 
بهینه‌یابی استفاده شد. در هر جزیره تعداد کل 
کروموزوم‌ها در هر نسل برابر با ۲۰ در نظر گرفته شده 
است لذا تعداد کل کروموزوم‌ها به‌ازای ۵۰۰ نسل برابر 
با ۱۰۰۰۰ می‌شود. تنها تحلیل یک کروموزوم بدون 
اعضاده از زاهکار رل لین سافتا طول نو کشسی ی 
به‌همین دلیل در تحقیقات گذشته کم‌تر به بهینه‌سازی 
چنین سازه‌هایی اقدام شده است. 

از لیست مقاطع زیر در امر بهینه‌یابی استفاده 
می‌شود: 


22, 24: 2.6, 28, 30, 32, 3.4, 6, 

_ [32.8 40, 42, 44, 46, 4 2 9 
5. 56 58 60 62 64 ۵ 
7.0, 7.2, 7.4, 76,78, 80, 82, 4 


تعداد کروموزوم‌هایی که شایستگی آن‌ها توسط 
عامل هو شمند فاز دوم مورد ارزیابی فسزا تین کر که 


سال بیست و چهارم شماره دو» ۱۳۹۲ 


عل ی کاوه- وحید رضا کلات جاری- پوریا منصوریان 


به‌طور متوسط برابر با ۵۵۰۰ است. بنابراین مشابه آن‌چه 
فو‌شال فی دک شتا مات یبن از نمی از 
کروموزوم‌ها توسط عامل هوشمند فاز دوم تعیین 
تقلیل می‌یابد. عدم تحلیل تکراری این کروموزوم‌ها 
همان گونه که مشاهده شل ثایر استفاده از روش 
جبری نیروها حدود ۹۵ درصد است و لذا در امر 
بهینه‌یابی از این روش به‌جای روش سختی استفاده 
شده است. هم‌چنین از همه‌ی عوامل هوشمند مصنوعی 
در این مثال استفاده شده اشنتت: 
ریبدت امه در نسل 
تحتانی احتمال خرابی سیستم سازه‌ای سازه‌ی بهینه 
به‌ترتیب براپر ۲ ۹/۸۰۷/۱۰۲ و 2/۷۹۷/۱۰۲ به‌دست 


حدول ۱۳ سطح مقطع سازه‌ی بهینه در مثال ۶ 


0۳/993 
00310 
سطح مقطع (سانتی‌متر مربع) 


۳/۲ ۱ 


شماره گروه 


1 
۳/۸ 
۳/۲ 
۳/۲ ۵ 
۳/۸ 1 
۳/۲ ۷ 


| - | حم 


وزن (18) 2۳۰۳۳ 


در این مثال برای تعیین زمان صرفه‌جویی شده در 
اثر استفاده از پسردازش موازی. از حالت ساده‌ی 
الگوریتم یعنی از یک پردازشگر استفاده شد و هر ۱۰۰ 
کروموزوم در یک جزیره بهینه شدند. زمان محاسبات 
در این حالت به کمی بیش‌تر از ۵ برابر حالت موازی 
افزایش یافت. دلیل این امر مشابه منال قبل, انجام 


سال بیست و چهارم. شماره دو. ۱۳۹۲ 


۳۳ 


سریع‌تر عملیاتی چون رتبه‌بندی کروموزوم‌ها؛ ادغام و 
برای ۲۰ کروموزوم نسبت به ۱۰۰ کروموزوم است و 
نتیجه‌ی به‌دست آهتله متشنلدا این است که زمان 
محاسبات در پردازش موازی جزیره‌ای رابطه‌ای 


معکوس و خطی با تعداد پردازشگرها دارد. 


نتیجه گیری 


در این تحقیق برای تسریع در روند تعیین مسیرهای 
خرابی محتمل و احتمال خرابی سیستم سازه ای سه 
راهکار و برای تسریع در امر بهینه‌یابی یک راهکار 
ارائه شد. 

در راهکار اول. ابتدا احتمال وقوع یسک مسیر 
خرابی با کوچک‌ترین احتمال مشترک هر یک از 
اعضای آن مسیر با اولین عضو تخمین زده شد. پس از 
آن از تولید شاخه‌های منشعب از عضوی که 
کوک تریه | ال متفرگ ترا هتعاس سیر گیرض 
گردید. دلیل اين امر وابستگی این مسیرها بوده‌است و 
تولید آن‌ها چیزی به اجتماع رویدادهای خرابی سازه 
اضافه نمی‌کند. جلوگیری از تولید این گونه مسیرها دو 
مزیت مهم به دنبال داشت. اول این‌که تعداد زیادی از 
مسیرهای همبسته‌ی دیگر در جمع احتمال مسیرهای 
خرابی شرکت نکرد و کران بالای کوچک‌تری برای 
احتمال خرابی سیستم سازه‌ای به‌دست آمد. مزیت دوم 
به‌طور بدیهی کاهش حجم و زمان محاسبات بود. 
استفاده از این راهکار علاوه بر باریک کردن کرانه‌های 
احتمال خرابی. زمان محاسبات را از ۵۰ تا نزدیک به 
۰ درصد کاهش داد. مجدداً بادآوری می‌شود که 
زمان تحلیل قابلیت اعتماد بدون استفاده از این راهکار 
معمولاً چندین ساعت است و لذا کاهش این زمان به 
چندین ثانیه به‌معنی کاهش حدود صد در صد است. 

در این تحقیق به‌عنوان راهکار دوم روشی سریع 
برای تحلیل خرپا با توپولوژی متغیر و استفاده از آن در 
تین قاپلیت اعتماه آزاند شل, این اسر با اتخاه ود 


نشریه مهندسی عمران فردوسی 


۳ 


در روابط روش جبری نیروها حاصل شد. تحلیل خرپا 
با روش مرسوم سختی نیازمند حل دستگاه معادلات با 
اندازه‌ی تعداد درجات آزادی فعال گرهی است در 
حالی که در روش نرمی اندازه‌ی دستگاه معادلات به 
فان قایمی ترا کاهشی مت بانن شر اکن رازه 
تعداد درجات نامعینی کوجک‌تر از تعداد درجات 
آزادی فعال گرهی است. مزیت دیگر روش جبری 
نیروها نسبت به روش مرسوم سختی این است که 
نیروها به طور مستقیم و بدون واسطه تعیین می‌شود. با 
توجه به دفعات بسیار زیاد تحلیل خرپا در تعیین 
مسیرهای خرابیء استفاده از روش نرمی پیشنهادی 
باعث صرفه جویی قابل توجه ۲۰ تا ٩۵‏ درصدی در 
زمان محاسبات شد. 

در این تحقیق به‌عنوان راهکار سوم از عامل 
هوشمند در دو فاز استفاده شد. از عامل هوشمند فاز 
اول در ذخیره‌سازی. شناسایی و کنترل مسیرهای 
خرابی استفاده شد که نتیجه‌ی آن جلوگیری از انجام 
هجو هو تک تون خامن سا سخاس ات 


حم 


بهینه یابی سریع حرپا با استفاده از الگوریتم ژنتیک ... 


بود. از عامل هوشمند فاز دوم برای شناسایی 
کروموزوم‌های شایسته‌ی تکراری و تعیین شایستگی 
آن‌ها استفاده شد. همان‌گونه که بیان شد تحلیل این 
گونه کروموزوم‌ها زمان بسیار بیش‌تری نسبت به سایر 
کروموزوم‌ها صرف می‌کند. مطابق نتایج شتسه 
تأثیر این عامل هوشمند از عامل فاز اول به‌مراتتب 
بیش‌تر بوده است. با استفاده از عامل هوشمند فاز دوم 
حدود ۷۰ درصد در زمان محاسبات صرفه‌جویی 
به‌دست آمد. 

با استفاده از راهکار چهارم ان تحقیق 
کروموزوم‌هایی که شایستگی کمی دارند در حین 
بل سدقم وایل, یل هی می سوبد: 
کروموزوم‌هایی که احتمال خرابی آن‌ها به مقدار قاببل 
توجهی با مقدار مجاز احتمال خرابی (قید طراحی) 
فاصله دارد در همان اوایل تحلیل شناخته می‌شوند و 
تحلیل آن‌ها زودتر از موعد مقرر پایان می‌یابد. استفاده 
از این راهکار باعث کوتاه شدن زمان تحلیل شد. 
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1 1 


(پاییز ۱۳۸۸). 
۲)0۲ 10 ,"۵006006و۳ ببن عصلونا فمصطاتنمعله ۵۶ وممتاه‌صا۳۵0 ویک رتاه‌وتصتدل رک رصماتامممع 
2004(۰) وکتمطوتامانظ )اما تعظ 20 دمم 
ماگ رصملاممتصصتامن میاه 00و02 واناصهنامز ملهمیمم نها و۱۷ رقتلق۳2۵۵202 .۷ رقلکلقصهم‌جطمو 1[ 
2004(۰) ,429-440 .00 و26 ومحتامن0 .ص01و1 ۱۷۲۵0 
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33. 


تشریه مهندسی عمران فردوسی سال بیست و چهارم شماره دو» ۱۳۹۲ 


